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Abstract: In a series of recently published papers we reanalyzed the exist- 
ing treatments of the Veneziano and Veneziano- like amplitudes and the models 
associated with these amplitudes. In this work we demonstrate that the already 
obtained new partition function for these amplitudes can be exactly mapped 
into that for the Polychronakos-Frahm (P-F) spin chain model which, in turn, 
is obtainable from the Richardon-Gaudin (R-G) XXX model. Reshetikhin and 
Varchenko demonstrated that such a model is obtainable as a leading approxi- 
mation in their WKB-type analysis of solutions of the Knizhnik-Zamolodchikov 
(K-Z) equations. The linear independence of solutions of these equations is con- 
trolled by determinants (discovered by Varchenko) whose explicit form up to 
a constant coincides with the Veneziano (or Veneziano-like) amplitudes. In 
the simplest case, when K-Z equations are reducible to the Gauss hypergeo- 
metric equation, the determinantal conditions coincide with those which were 
discovered by Kummer in 19-th century. Kummer's results admit physical in- 
terpretation crucial for providing needed justification associating determinantal 
formula(s) with Veneziano-like amplitudes. General results arc illustrated by 
many examples. These include but are not limited to only high energy physics 
since all high energy physics scattering processes can be looked upon from much 
broader stochastic theory of random fragmentation and coagulation processes 
recently undergoing active development in view of its applications in disciplines 
ranging from ordering in spin glasses and population genetics to computer sci- 
ence, linguistics and economics, etc. In this theory Veneziano amplitudes play 
a central (universal) role since they are the Poisson-Dirichlet-type distributions 
for these processes (analogous to the more familiar Maxwell distribution for 
gases). 

Keywords: Polychronakos and Richardson-Gaudin spin chains, Knizhnik- 
Zamolodchikov equations, determinantal formulas, Veneziano amplitudes, ran- 
dom fragmentation-coagulation processes. 
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1 Introduction 



Since time when quantum mechanics (QM) was born (in 1925-1926) two seem- 
ingly opposite approaches for description of atomic and subatomic physics were 
proposed respectively by Heisenberg and Schrodinger. Heisenberg's approach is 
aimed at providing an affirmative answer to the following question: Is combina- 
torics of spectra (of obsevables) provides sufficient information about micro- 
scopic system so that dynamics of such a system can be described in terms of 
known macroscopic concepts? Schrodinger's approach is exactly opposite and is 
aimed at providing an affirmative answer to the following question: Using some 
plausible mathematical arguments is it possible to find an equation which under 
some prescribed restrictions will reproduce the spectra of observables? Although 
it is widely believed that both approaches are equivalent, already Dirac in his 
lectures on quantum field theory [1] noticed (without much elaboration) that 
Schrodinger's description of QM contains a lot of "dead wood" which can be 
safely disposed altogether. According to Dirac "Heisenberg's picture of QM is 
good because Heisenberg's equations of motion make sense" . 

To our knowledge, Dirac's comments were completely ignored, perhaps, be- 
cause he had not provided enough evidence making Heisenberg's description of 
QM superior to that of Schrodinger's. In recent papers [2,3] we found examples 
supporting Dirac's claims. From the point of view of combinatorics, there is 
not much difference in description of QM, quantum field theory and string the- 
ory. Therefore, in this paper we choose the Heisenberg's point of view on string 
theory using results of our recent works in which we re analyzed the existing 
treatments connecting Veneziano (and Vcneziano-like) amplitudes with the re- 
spective string-theoretic models. As result, we were able to find new tachyon-free 
models reproducing Veneziano (and Veneziano-like) amplitudes. In this work 
the result of our papers [4-6] which will be called as Part I, Part II and Part 
III respectively, are developed further to bring them in correspondence with 
those proposed by other authors. Without any changes in the already devel- 
oped formalism, we were able to connect our results with an impressive number 
of string-theoretic models, including the most recent ones. Nevertheless, be- 
low we argue that, although physically plausible, the established connections 
(in the way they are typically treated in physics literature) are mathemati- 
cally ill founded. To correct this deficiency, in Section 5 we use some works by 
mathematicians. Of particular importance for us are the works by Reshetikhin 
and Varchenko [7] and by Varchcnko summarized in Varchcnko's MIT lec- 
ture notes [8] . These works enabled us to relate Veneziano (and Veneziano-like) 
amplitudes (e.g. those describing tttt scattering) to Knizhnik-Zamolodchikov 
(K-Z) equations and, hence, to WZNW models. This is achieved by employing 
known connections between the WZNW models and spin chains. In the present 
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case, between the K-Z equations and the XXX-type Richardson-Gaudin mag- 
netic chains as described in Section 5. Sections 2-4 contain mathematically less 
sophisticated results aimed at providing needed physical motivations and back- 
ground. For this purpose in section 2 we replaced mathematically sophisticated 
derivation of the Veneziano partition function by considerably simpler combi- 
natorial derivation of such function. As a by product of this effort we were able 
to uncover the connections with spin chains already at this stage of our inves- 
tigation. To strengthen this connection, in Section 3 we demonstrate that the 
obtained Veneziano partition function coincides with the Polychronakos-Frahm 
(P-F) partition function for the ferromagnetic spin chain model. Although such 
a spin chain was extensively studied in literature, we discuss different paths in 
Section 4 aimed at establishing links between the P-F spin chain and variety 
of string-theoretic models, including the most recent ones. This is achieved 
by mapping combinatorial and analytical properties of the P-F spin chains 
into analogous properties of spin chains used for description of the stochas- 
tic process known as asymptotic simple exclusion process (ASEP). To make 
our presentation self-contained, we provide in Appendix A basic information 
on ASEP sufficient for understanding the results discussed in the main text. 
In addition, in the main text we provide some information on Kardar-Parisi- 
Zhang (KPZ) and Edwards- Wilkinson (EW) equations which are just different 
well defined macroscopic limits of the microscopic ASEP equations. We do this 
with purpose of reproducing variety of string-theoretic models, including the 
most recent ones. Such a success have not deterred us from looking at other, 
more rigorous (mathematically) approaches. These are discussed in Sections 5 
and, in part, in Section 6. These sections are interrelated and contain the most 
important results of this paper. While the content of Section 5 was already 
briefly discussed, the content of Section 6 provides the strongest independent 
support to the results and conclusions of Section 5. At the same time, this 
section can be read independently of the rest of the paper since it contains 
some important facts from the theory of random fragmentation and coagu- 
lation processes [9-11] which is currently in the process of rapid development 
because of its wide applications ranging from theory of spin glasses and pop- 
ulation genetics to computer science, linguistics and economics, etc. In high 
energy physics this theory was developed for some time by Mekjian, e.g. see 
[12] and references therein. Since our Section 6 is not a review, our treatment 
of topics discussed in it is markedly different from that developed in Mekjian's 
papers and is subordinated to the content of Section 5. Specifically, the main 
result of Section 5 is the deteminantal formula, equation (5.47), which up to a 
constant coincides with the Veneziano (or Veneziano-like) amplitude. A special 
case of this formula produces known pion-pion scattering amplitude. In Section 
6 we argue that: 1. Veneziano amplitudes play the central role in the theory of 
random fragmentation and coagulation processes where they are known as the 
Poisson-Diriclct (P-D) probability distributions. 2. The discrete spectra of all 
exactly solvable quantum mechanical (QM) problems can be rederived in terms 
of some P-D stochastic processes. This is so because all exactly solvable QM 
problems involve some kind of orthogonal polynomials-all derivable from the 
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Gauss hypergcomctric function which admits an interpretation in terms of the 
P-D process. 3. Since in the simplest case the K-Z equations are reducible to 
the hypergeometric equations, the processes they describe are also of P-D type. 
4. In the case of Gauss hyprgeometric equation, the determinantal formula 
(5.47) is reduced to that obtained by Kummer in 19th century. To facilitate 
understanding and appreciation of these facts and to demonstrate utility of 
the obtained results beyond the scope of high energy physics, in Section 6 we 
discuss some applications of the developed formalism to genetics and chemical 
kinetics. 

2 Combinatorics of Veneziano amplitudes and 
spin chains: qualitative considerations 

In Part I, we noticed that the Veneziano condition for the 4-particle amplitude 
given by 

a(s) +a(t) + a(u) = -1, (2.1) 

where a(s), a(t), a(u) € Z, can be rewritten in more mathematically suggestive 
form. To this purpose, following [13], we need to consider additional homogenous 
equation of the type 

a{s)m + a(t)n + a(u)l + k ■ 1 = (2.2) 
with m, n, l 7 k being some integers. By adding this equation to (2.1) we obtain, 

a{s)fh + a(t)n + a(u)l = k (2.3a) 

or, equivalently, as 

m + n 2 + n 3 = N, (2.3b) 

where all entries by design are nonnegative integers. For the multiparticle case 
this equation should be replaced by 

n + --- + n k =N (2.4) 

so that combinatorially the task lies in finding all nonnegative integer combina- 
tions of no, ...,n,h producing (2.4). It should be noted that such a task makes 
sense as long as N is assigned. But the actual value of N is not fixed and, hence, 
can be chosen quite arbitrarily. Equation (2.1) is a simple statement about the 
energy -momentum conservation. Although the numerical entries in this equa- 
tion can be changed as we just explained, the actual physical values can be 
subsequently re obtained by the appropriate coordinate shift. Such a procedure 
should be applied to the amplitudes of conformal field theories (CFT) with 
some caution since the periodic ( or antiperiodic, etc.) boundary conditions 
cause energy and momenta to become a quasi -energy and a quasi momenta (as 
it is known from solid state physics). 
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The arbitrariness of selecting N reflects a kind of gauge freedom. As in other 
gauge theories, we may try to fix the gauge by using some physical consider- 
ations. These include, for example, an observation made in Part I that the 4 
particle amplitude is zero if any two entries into (2.1) are the same. This fact 
prompts us to arrange the entries in (2.3b) in accordance with their magnitude, 
i.e. m > n 2 >ri3. More generally, we can write: n > n% > ■ ■ ■ > nu > 10. 

In Section 6 we demonstrate that if the entries in this sequence of inequalities 
arc treated as random nonnegative numbers subject to the constraint (2.4), these 
constrains are necessary and sufficient for recovery of the probability density for 
such set of random numbers. This density is known in mathematics as Dirichlet 
distribution^ [9-11,14]. Without normalization, integrals over this distribution 
coincide with Veneziano amplitudes. 

Provided that (2.4) holds, we shall call such a sequence a partition and shall 
denote it as n =(no, ...,rik). If n is partition of N, then we shall write n h N. 
It is well known [15,16] that there is one- to -one correspondence between the 
Young diagrams and partitions. We would like to use this fact in order to design 
a partition function capable of reproducing the Veneziano (and Veneziano-like) 
amplitudes. Clearly, such a partition function should also make physical sense. 
Hence, we would like to provide some qualitative arguments aimed at convincing 
our readers that such a partition function does exist and is physically sensible. 

We begin with observation that there is one- to- one correspondence between 
the Young tableaux and directed random walkfl It is useful to recall details of 
this correspondence now. To this purpose we need to consider a square lattice 
and to place on it the Young diagram associated with some particular partition. 

Let us choose some n x rh rectangl^ so that the Young diagram occupies 
the left part of this rectangle. We choose the upper left vertex of the rectangle 
as the origin of the xy coordinate system whose y axis (South direction) is 
directed downwards and x axis is directed Eastwards. Then, the South-East 
boundary of the Young diagram can be interpreted as directed (that is without 
self intersections) random walk which begins at (0, — in) and ends at (n,0). 
Evidently, such a walk completely determines the diagram. The walk can be 
described by a sequence of 0's and l's. Say, for the x— step move and 1 for the 
y— step move. The totality Af of Young diagrams which can be placed into such 
a rectangle is in one-to-one correspondence with the number of arrangements 
of 0's and l's whose total number is rh + n. Recalling the Fermi statistics, the 
number Af can be easily calculated and is given by Af = (m + n)!/m!n|E It can 

1 The last inequality: nj. > 1, is chosen only for the sake of comparison with the existing 
literature conventions, e.g. see Ref. [15]. 

2 For reasons explained in Section 6 it is also called the Poisson-Dirichlet distribution. 

3 Furthermore, it is possible to map bijectively such type of random walk back into Young 
diagram with only two rows, e.g. read [17], page 5. This allows us to make a connection with 
spin chains at once. In this work we are not going to use this route to spin chains in view of 
the simplicity of alternative approaches discussed in this section. 

4 Parameters h and rh will be specified shortly below. 

5 We have suppressed the tildas for n and m in this expression since these parameters are 
going to be redefined below anyway. 
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be represented in two equivalent ways 

(n + l)(n + 2) • • • (n + to) _ 



(to + n)\/m\n\ 



(to + l)(m + 2) • • • (n + to) _ 



n + to 
m 

to + n 
n 



(2.5) 



Let now p(iV; k, to) be the number of partitions of N into < k nonnegative 
parts, each not larger than to. Consider the generating function of the following 
type 



F(k,m\q)=}^p(N;k,m)q N (2.6) 

JV=0 

where the upper limit S will be determined shortly below. It is shown in 

k + m _ k + to 

TO k 



Refs. [15,16] that T{k, to | g) = 



where, for instance. 



k + to 
m 



E| From this result it should be clear that the expression 

k + n 



k + to 

TO 



+ TO 

rn 



9=1 



is the g— analog of the binomial coefficient I '" ' "" J . In literature [15,16] this 
q— analog is known as the Gaussian coefficient. Explicitly, it is defined as 

(q a - l){q a ~ l -!)••• (q a - b+1 - 1) 



b-l 



!)•• •(<?-!) 



(2.7) 



for some nonegative integers o and b. From this definition we anticipate that 
the sum defining generating function !F{k, m \ q) in (2.6) should have only finite 
number of terms. Equation (2.7) allows easy determination of the upper limit S 
in the sum (2.6). It is given by km. This is just the area of the k x to rectangle. 
In view of the definition of p(N; k, m), the number to = N — k. Using this fact 

r N I 

(2.6) can be rewritten as: F{N, k | q) = , .This expression happens to be 

L J q 

the Poincare' polynomial for the Grassmannian Gr(m, k) of the complex vector 
space C^of dimension N as can be seen from page 292 of the book by Bott and 
Tu, [180- From this (topological) point of view the numerical coefficients, i.e. 
p{N;k,m), in the q expansion of (2.6) should be interpreted as Betti numbers 
of this Grassmannian. They can be determined recursively using the following 



6 On page 15 of the book by Stanley [16], one can find that the number of solutions N(n, k) 

( Ti 4- k 1 

in positive integers to yi + ... + = n + k is given by 



solutions in nonnegative integers to ii 



I , while the number of 
n + k \ 

- it = n is | , J . Careful reading of Page 

15 indicates however that the last number refers to solution in nonnegative integers of the 
equation xq + ... + xj, = n. This fact was used essentially in (1.21) of Part I. 

7 To make a comparison it is sufficient to replace parameters t 2 and n in Bott and Tu book 
by q and N. 
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property of the Gaussian coefficients [4], page 26, 



' n + 1 ' 




n 




n 


k + l 


g 


k + l 


+ q n-k 

q 


k 



and taking into account that 



n 




(2- 



= 1. We refer our readers to Part II for 



mathematical proof that T{N, k \ q) is indeed the Poincare' polynomial for the 
complex Grassmannian. With this fact proven, we notice that, due to relation 
m = N — k, it is sometimes more convenient for us to use the parameters m 
and k rather than N and k. With such a replacement we obtain: 



T(k, m 



k 

n 1 ! 



k + m 
k 

n m-\-i 



(q 



!)(<? 



k-\-m— 1 



m+1 



i) 



(q k -l)(q 



k-1 



l)---(g-l) 



(2.9) 



This result is of central importance. In our work, Part II, considerably more 
sophisticated mathematical apparatus was used to obtain it (e.g. see equation 
(6.10) of this reference and arguments leading to it). 

In the limit : q — > 1 (2.9) reduces to N as required. To make connections 
with results known in physics literature we need to re scale q's in (2.9), e.g. 
let q = t~ . Substitution of such an expression back into (2.9) and taking the 
limit t — > 1 again produces Af in view of (2.5). This time, however, we can 
accomplish more. By noticing that in (2.4) the actual value of TV deliberately 
is not yet fixed and taking into account that m = N — k we can fix N by fixing 
m. Specifically, we would like to choose m = 1 • 2 • 3 • • • k and with such a choice 
we would like to consider a particular term in the product (2.9), e.g. 



S(i) 



1-t 



i+- 



l-t 



(2.10) 



In view of our "gauge fixing" the ratio m/i is a positive integer by design. This 
means that we are having a geometric progression. Indeed, if we rescale t again 
: t — > t 2 , we then obtain: 



S(i) 



l + t 2 + ---+t 



2m 



(2.11) 



with m = y . Written in such a form the above sum is just the Poincare' poly- 
nomial for the complex projective space CP m . This can be seen by comparing 
pages 177 and 269 of the book by Bott and Tu [18]. Hence, at least for some 
m's, the Poincare' polynomial for the Grassmannian in just the product of the 
Poincare' polynomials for the complex projective spaces of known dimensional- 
ities. For m just chosen, in the limit t — > 1, we reobtain back the number Af as 
required. This physically motivating process of gauge fixing just described can 
be replaced by more rigorous mathematical arguments. The recursion relation 
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(2.8) introduced earlier indicates that this is possible. The mathematical de- 
tails leading to factorization which we just described can be found, for instance, 
in the Ch-3 of lecture notes by Schwartz [19]. The relevant physics emerges by 
noticing that the partition function Z(J) for the particle with spin J is given 
by [20] 

Z(J) = tr{e~ pH ^) = e cJ + e c(J -^ + ■■■ + e~ cJ 

= e cJ {l + e- c + e- 2c + --- + e- 2cJ ), (2.12) 

where c is known constant. Evidently, up to a constant, Z(J) ~ S(i). Since 
mathematically the result (2.12) is the Weyl character formula, this fact brings 
the classical group theory into our discussion. More importantly, because the 
partition function for the particle with spin J can be written in the language of 
N=2 supersymmetric quantum mechanical models, as demonstrated by Stone 
[20] and others [21], the connection between the supersymmetry and the classical 
group theory is evident. It was developed in Part III. 

In view of arguments presented above, the Poincare' polynomial for the 
Grassmannian can be interpreted as a partition function for some kind of a spin 
chain made of apparently independent spins of various magnitude^. These 
qualitative arguments we would like to make more mathematically and physi- 
cally rigorous. The first step towards this goal is made in the next section. 

3 Connection with the Polychronakos-Frahm spin 
chain model 



The Polychronakos-Frahm (P-F) spin chain model was originally proposed by 
Polychronakos and described in detail in [23]. Frahm [24] motivated by the 
results of Polychronakos made additional progress in elucidating the spectrum 
and thermodynamic properties of this model so that it had become known as 
the P-F model. Subsequently, many other researchers have contributed to our 
understanding of this exactly integrable spin chain model. Since this paper is 
not a review, we shall quote only works on P-F model which are of immediate 
relevance. 

Following [23] , we begin with some description of the P-F model. Let <jf (a — 
1, 2, n 2 — 1) be SU(n) spin operator of i-th particle and let the operator <7j.,- 
be responsible for a spin exchange between particles i and j, i.e. 

<*i = ^ (3-1) 



8 We hope that no confusion is made about the meaning of N in the present case. 

9 In such a context it can be vaguely considered as a variation on the theme of the Polyakov 
rigid string (Grassmann a model, Ref. [22], pages 283-287), except that now it is exactly 
solvable in the qualitative context just described and, below, in mathematically rigorous 
context. 
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In terms of these definitions, the Calogero-type model Hamiltonian can be writ- 
ten as [25,26] 

n = lz(pl + ^) + z ! (l ~ aij l (3-2) 

1 i i<3 [Xi - Xj) 

where I is some parameter. The P-F model is obtained from the above model 
in the limit I — > ±oo . Upon proper rescaling of TL in (3.2), in this limit one 
obtains 

Hv-r = -sign{l) £ aij 2 , (3.3) 

i<j \ x i ~ Xj) 

where the coordinate Xi minimizes the potential for the rescaled Calogero modef^l. 
that is 

i<j (Xi - Xj) 

It should be noted that Ti-p^jr is well defined without such a minimization, that 
is for arbitrary real parameters Xi . This fact will be further explained in Section 
5. In the large I limit the spectrum of Ji is decomposable as 

E = E C + lEp-r (3.5) 

where Eq is the spectrum of the spinless Calogero model while E-p^jr is the 
spectrum of the P-F model. In view of such a decomposition, the partition 
function for the Hamiltonian 7i at temperature T can be written as a product: 
Z H {T) =Z c (T)Z-p^jr(T/l). From here, one formally obtains the result: 

Zr-AT) = lim §4^?- ( 3 - 6 ) 

It implies that the spectrum of the P-F spin chain can be obtained if both the 
total and the Calogero partition functions can be calculated. In [23] Polychron- 
akos argued that Zq(T) is essentially a partition function of N noninteracting 
harmonic oscillators. Thus, we obtain 

N 1 i 

Z C (N;T) = n 1 q = exp(-/3),/3= {k B Ty l . (3.7) 

i=l 1 — 1 

Furthermore, the partition function Zu (T) according to Polychronakos can be 
obtained using Zc(N;T) as follows. Consider the grand partition function of 
the type 

oo / oo \ n 

5 = £ Z n (N;T)y N = £ Z c (L;T)y L (3.8) 

N=0 \L=Q / 

where n is the number of flavorJ^l Using this definition we obtain 

n 

Z n {N;T)= J2 X[Zc(h:T). (3.9) 



10 The Calogero model is obtainable from the Hamiltonian (3.2) if one replaces the spin 
exchange operator Uij by 1. Since we are interested in the large I limit, one can replace the 
factor 1(1 — 1) by I 2 in the interaction term. 

llr That is n the same number as n in SU(n). 
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Next, Polychronakos identifies Z n (N;T) with Z n (T). Then, with help of (3.6) 
the partition function Z-p-jr(T) is obtained straightforwardly as 

ft(W) 

Z V _?{N-T)= ]T ^4 . (3.10) 

ri(i-? r ) 

i=l r=l 

Consider this result for a special case: n = 2. It is convenient to evaluate the 
ratio first before calculating the sum. Thus, we obtain: 

N 

fi (W) _ (i-g)---(i-O -n"H ^ 

n n(i-r) 

2=1 r=l 

(3.H) 

where the Poincare' polynomial T{N, k \ q) for the Grassmanian of the complex 
vector space C" of dimension N was obtained in the previous section. Indeed 
(3.11) can be trivially brought into the same form as given in our equation (2.9) 
using the relation m + k = N . To bring (2.9) in correspondence with equation 
(4.1) of Polychronakos [23], we use the second equality (2.9) in which we make 
a substitution: m = N — k. After this replacement, (3.10) acquires the following 
form 

N k 1 - /-•+! 
fe=0 i=0 1 

coinciding with equation (4.1) by Polychronakos. This equation corresponds to 
the ferromagnetic version of the P-F spin chain model. To obtain the antifcrro- 
magnetic version of the model requires us only to replace q by q -1 in (3.12) and 
to multiply the whole r.h.s. by some known power of q. Since this factor will 
not affect thermodynamics, following Frahm [24], we shall ignore it. As result, 
we obtain 

AT k 1 N—i-\-l 

Z«l A N- D = £ <Z W2 - fc)2 II \ q _ ai . (3-13) 

k=0 i=0 ^ 

in accord with Frahm's equation (21). This result is analyzed further in the 
next section. 
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4 Connections with WZNW model and XXX 
s=l/2 Heisenberg antiferromagnetic spin chain 

4.1 General remarks 



To establish these connections we follow work by Hikami [27]. For this purpose, 
we introduce the notation 



i=l 

allowing us to rewrite (3.13) in the equivalent form 

N k <v_, i i N 



1 ~Q N 1+1 _ V- JN/2-k) 2 (g)jy 



Z?_ A N; T) = £ II = E 



( 1 



. (4-2) 

Consider now the limiting case (N — > oo) of the obtained expression. For this 
purpose we need to take into account that [andrewsl] 



lim 

N^oo 



N 
k 



To use this asymptotic result in (4.2) it is convenient to consider separately 
the cases of N being even and odd. For instance, if N is even, we can write: 
N = 2m. In such a case we can introduce new summation variables: I = k — m 
and/or I = m — k. Then, in the limit TV — > oo (that is m^ oo) we obtain 
asymptotically 

Z a J_r{oo;T) = — /• (4-4a) 

\Q)oO i = -oo 

in accord with [27]. Analogously, if N = 2m + 1, we obtain instead 



Z. 



af 



-. oo 

(oo; T) = — Y, 1 {l+h)2 - (4.4b) 



According to Melzer [28] and Kedem, McCoy and Melzer [29], the obtained 
partition functions coincide with the Virasoro characters for SUi(2) WZNW 
model describing the conformal limit of the XXX (s=l/2) antiferromagnetic 
spin chain [30]. Even though equations (4.4a) and (4.4b) provide the final 
result, they do not reveal their physical content. This task was accomplished 
in part in the same papers where connection with the excitation spectrum of 
the XXX antiferromagnetic chain was made. To avoid repetitions, below we 
arrive at these results using different arguments. By doing so many new and 
unexpected connections with other stochastic models will be uncovered. 
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4.2 Method of generating functions and q-deformed har- 
monic oscillator 

We begin with definitions. In view of (2. 9), (3. 12) and (4.2), we would like to 
introduce the Galois number Gn via 



N 



fc=0 



N 
k 



(4.5) 



This number can be calculated recursively as it was shown by Goldman and 
Rota [31] with the result 



Gjv+i — 2G 



JV 



1 N 1) Gjv-i- 



(4.6) 



Alternative proof is given by Kac and Cheung [32]. To calculate Gat we have 
to take into account that Go = 1 and Gi = 2. These results can be used 
as a reference when one attempts to calculate the related Rogers-Szego (R-S) 
polynomial H^(t) defined as [33] 



JV 



H N {t;q) :=Y, 



fe=0 



N 
k 



(4.7) 



so that Hn(1) — G/vO Using [32] once again, we find that ifjv(i) obeys the 
following recursion relation 



H N+1 (t) = (l+t)H N (t) + (q N - 1) tH N -x(t) 



(4.8) 



which for t = 1 coincides with (4.6) as required. The above recursion relation 
is supplemented with initial conditions. These are : Hq — 1 and Hi = 1 + 1. 

At this point we would like to remind our readers that for t — 1 according 
to (3.12) we obtain: Zj,_jr(N ;T) = Gn- Hence, by calculating Hpj(t) we shall 
obtain the partition function for the P-F chain. 

To proceed with such calculations, we follow [34]. In particular, consider 
first an auxiliary recursion relation for the Hermite polynomials: 



H n+ i(x) = 2xH n (x) - 2nH n ^ 1 {x) 
supplemented by the differential relation 



d 

dx 



H n (x) = 2nH n -i{x) 



(4.9a) 



(4.9b) 



which, in view of (4.9a), can be conveniently rewritten as 



H n+ x(x) = (2x - 



dx 



)H n (x). 



(4.9c) 



2 For brevity, unless needed explicitly, we shall suppress the argument q in Hf^(t;q). 
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This observation prompts us to introduce the raising operator R = 2x — so 

dx 

that we obtain: 

R n H Q {x)=H n (x). (4.10) 
The lowering operator can be now easily obtained again using (4.9). We get 
1 d 



2 dx 



H n (x) = LH n (x) = nH n -i{x) 



(4.11) 



so that [L, R] = 1 as required. Based on this,the number operator N can be 
obtained as N = RL so that NH n (x) = nH n (x) or, explicitly, using provided 
definitions, we obtain 



(^-2xl + 2n)H n{ x)=0. 



n + 1 >, L I n >- 



(4.12) 
1 > and , 



Evidently, we can write: R \ n > 
< m | n >= n\8 mn as usual. 

We would like now to transfer all these results to our main object of interest 
the recursion relation (4.8). To this purpose, we introduce the difference oper 
ator A via 

AH N (t) :=H N (t)-H N (qt). 
Using definition (4.7) we obtain now 



AH N (t) = (l-q N )tH N _ 1 (t) 
where we took into account that 



(4.13) 
(4.14) 



" N ' 




N 


k 


i 


N-k 



Using this result in (4.8) we obtain at once 

H N+1 (t) = [(l+t)-A]H N (t). 



(4.16) 



(4.17) 



This, again, can be looked upon as a definition of a raising operator so that we 
can formally rewrite (4.17) as 



nH N (t) = H N+1 (t). 

The lowering operator can be defined now as 



C := -A 

x 



so that 



(4.18) 

(4.19) 
(4.20) 



CH N (t) = {l-q N )H N _ 1 {t). 
The action of the number operator Af = 1ZC is now straightforward, i.e. 

AfH N (t) = (l-q N )H N (t). (4.21) 
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Following Kac and Cheung [32] we introduce the q— derivative via 
By combining this result with (4.13) we obtain, 

D - f ^ = W=TY <4 ' 23 » 

This allows us to rewrite the raising and lowering operators in terms of q— derivatives. 
Specifically, we obtain: 

TZ:=(l + t)-{l-q)tD q (4.24) 

and 

I := D q . (4.25) 

While for the raising operator rewritten in such a way equation (4.18) still holds, 
for the lowering operator C we now obtain: 

CH N {t) = IzjF-HN-^t) = [N\H N -i{t). (4.26) 

The number operator N q is acting in this case as 

N q H N {t) = [N]H N (t). (4.27) 

We would like to connect these results with those available in literature on 
q— deformed harmonic oscillator. Following Chaichan et al [35], we notice that 
the undeformed oscillator algebra is given in terms of the following commutation 
relations 

aa + - a + a = 1 (4.28a) 
[N, a] = -a (4.28b) 

and 

[N,a+}=a+. (4.28c) 

In these relations it is not assumed a priori that N = a + a and, therefore, 
this algebra is formally different from the traditionally used [a, a + ] = 1 for 
the harmonic oscillator. This observation allows us to introduce the central 
element Z — N — a + a which is zero for the standard oscillator algebra. The 
deformed oscillator algebra can be obtained now using equations (4.28) in which 
one should replace (4.28a) by [floreanni&vinet] 



aa 



+ -qa + a=l. (4.28d) 



Consider now the combination K := ClZ-qAf acting on Hn using previously 
introduced definitions. A simple calculation produces an operator identity: 
CR — qj\f =1 so that we can formally make a provisional identification : C -^a 
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and 1Z -^a + . To proceed, we need to demonstrate that with such an identifica- 
tion equations (4.28 b,c) hold as well. For this to happen, we should properly 
normalize our wave function in accord with known procedure for the harmonic 

oscillator where we have to use | n >= —= {a + ) n | >. In the present case, we 

Vn! 

have to use I N >= — ■ (lZ) n I > as the basis wavcfunction while making 

an identification: | >= H (t). The eigenvalue equation (4.27), when written 
explicitly, acquires the following form: 

" + = «• ( 4 - 29 ) 

4.3 The limit q — > 1^ and emergence of the Stieltjes-Wigert 
polynomials 



Obtained results need further refinements for the following reasons. Although 
the recursion relations (4.8), (4.9) look similar, in the limit q — > l* (4.8) is 
not transformed into (4.9). Accordingly, (4.29) is not converted into equation 
for the Hermitc polynomials known for harmonic oscillator. Fortunately, the 
situation can be repaired in view of recent paper by Karabulut [37] who spotted 
and corrected some error in the influential earlier paper by Macfarlane [39]. 
Following [38] we define the translation operator T(s) as T(s) :— e s a^. Using 
this definition, the creation a^ 1 and annihilation a operators are defined as follows 

a f = * [q x+ i - {s)]T^ (s) (4.30a) 

v l _ 1 

where T^(s) = and, accordingly, 

a = / T?(s)[g x+ i - Ti{s)\. (4.30b) 
v 1 — 1 

Under such conditions, the inner product is defined in the standard way, that is 



DC 



(/,<?)= / r(x)g(x)dx (4.31) 



so that (q x y = q x and (d/dxY = —(d/dx) thus making the operator to be 
a conjugate of a in a usual way. The creation-annihilation operators just defined 
satisfy commutation relation (4.28 d). At the same time, the combination a^a 
while acting on the wave functions ^ n (to be defined below) produces equation 
similar to (4.27), that is 

a f a*„ = [n]* n = A„*„ . (4.32) 
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Furthermore, it can be shown that 

a*n = \An*n-i and a f *„ = y/X n +i^ n +i 



(4.33) 



in accord with previously obtained results. Next, we would like to obtain the 
wave function \t„ explicitly. To this purpose we start with the ground state 
a* = and use (4.30b) to get (for s=l/2j3 the following result 



(4.34) 



Let w(x) be some yet unknown function. Then, it is appropriate to look for 
solution of (4.34) in the form 



i x ) = const ■ w(x)q x 



(4.35a) 



provided that the function w[x) is periodic: w(x) — w(x+l/2). The normalized 
ground state function acquires then the following look 



* (x) = a w w(x)q x 
where the constant a w is given by 



(4.35b) 



/ oo 






I/* 


q x w(x) 




\— oo 







(4.35c) 



Using this result, can be constructed in a standard way through use of the 
raising operators. There is, however, a faster way to obtain the desired result. 
To this purpose, in view of (4.35b), suppose that ^„ (x) can be decomposed as 
follows 

a w w(x) 



*n (x) 



J2C'k(q)(-l) k q (n - k)/2 q {x - 



(4.36) 



fc=0 



where (q, q) n = (1 — q)(l — q 2 ) ■ ■ ■ (1 — q n ) and C^{q) is to be determined as 
follows. By applying the operator a) / U A n +i to (4.36) and taking into account 
that T~5 w(x) = w(x — 1/2) = w(x) (in view of the periodicity of w(x)) we 
end up with the recursion relation for C£(q): 



C- +1 (q)^q"Cj:(q) + Cj:_ 1 (q)- 



(4.37) 



This relation should be compared with that given by (2.8). Andrews [33], page 
35, demonstrated that (2.8) and (4.37) are equivalent. Hence, we obtain: 



n 
k 



(4.38) 



13 The rationale for choosing s=l/2 is explained in the same reference. 



17 



This implies that, indeed, up to a constant, the obtained wavefunction should 
be related to the Rogers-Szego polynomial. This relation is nontrivial however. 
We would like to discuss it in some detail now. 

Following [37,39], let q — e~ c , where c is some nonegative number. Intro- 
duce the distributed Gaussian polynomials via 

oo 

* n (x) = ^(q)(-l) V^V*^ 2 . (4.39) 

fc=0 

These polynomials satisfy the following orthogonality relation: 



$ n (x)$ m (x)dx = ||*„(a;)|| 2 S mn (4.40) 



withal 

\\®n(x)\\ = (^y q -%V(wU- (4-41) 

This result calls for change in normalization of $„(a;), i.e.,<f> n (x) = j^rj^- 
Under such conditions <fi n (x) coincides with ^> n (x), provided that w(x) = 1. 
Introduce new variable: u = q~ 2x , and consider a shift: $ ra (x) — > 3> n (x — s). 
Using (4.39), we can write 

$ n (x-s) = u s exp{- (\nuf / (-4:\nq)}P n {u; s) (4.42a) 

where 

oo 

P n (u; s)=J2 CJt(q)(~l) k q- k / 2 q^ 2 u k . (4.42b) 
k=0 

The orthogonality relation (4.40) is converted then into 



duu^expi-ilnu) 2 /(-4\nq)}P n (u;s)P m (u;s) = 6 mn . (4.43) 



In view of (4.43) consider now a special case: s = 1/2. Then, the weight func- 
tion is known as lognormal distribution and polynomials P n {u; 1/2) (up to a 
constant ) are known as Stieltjes-Wigert (S-W) polynomials. Their physical rel- 
evance will be discussed below in Subsection 4.6. In the meantime, we introduce 
the Fourier transform of f(x) in the usual way as 



dxexp{2Tri6x)f(x) = f{6) {AAA) 



14 Notice that a = ( f (teg 2 * 2 J = (^r) 
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Then, the Parseval relation implies: 

oo oo 

[ $ n (x)$ m (x)dx = J <Z> n (O)$ m (0)dx = \\$ n (x)\\ 2 6 mn , (4.45a) 



causing 



l oo 

^(0) = (J) 3 exp(-( 7 r/ C )0 2 )^C?(g)(- (Z -^ 2 " 9 ) fc . (4.45b) 



fe=0 



Comparison between these results and (4.7) produces 

oo 

J H n (- q -h- 2 ™ e ; q )H m (- q -h- 2 ™ e ;q)cM-2(*/c)0 2 ) = (£) ' q~ n (q, q) n 6 r , 

— CO 

(4.46a) 

which can be alternatively rewritten as 



l 

j I ! "< UM )U,„[ -q 4, -'") l ),(2n0:q)<10 = q~"(q.q)„6 (4.401.) 



DC 



with $ 3 (6, q) = X) 1 n /2e • That is ^3 is one of the Jacobi's theta tunc- 

n— — 00 

tions. In order to use the obtained results, it is useful to compare them against 
those, known in literature already, e.g. see [40]. Equation (4.46a) is in agree- 
ment with (5) of [40] if we make identifications: k = n and c = \/2k, where k is 
the parameter introduced in this reference. With help of such an identification 
we can proceed with comparison. For this purpose, following 40] we introduce 
yet another generating function 



S n (t;q) :=Y, 



n 
k 



q k t k (4.47) 

k=o L ^ -I q 

so that the S-W polynomials can be written now as [41], page 197, 

S n (t;q) = (-ljVCv^W -1 ^*;^) = (-I)"?" (V^qUr'S^-q^q), 

(4.48) 

provided that < q < 1. Comparison between generating functions (4.7) and 
(4.47) allows us to write as well 

S n {t; q- 1 ) = H n (tq- n ;q), or, equivalently, H n (t; q' 1 ) = S n (q- n t: q) (4.49) 

Using this result we can rewrite the recursion relation (4.8) for H n (t; q) in terms 
of the recursion relation for S n (t;q) if needed and then to repeat all the argu- 
ments with creation and annihilation operators, etc. For the sake of space, we 
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leave this option as an exercise for our readers. Instead, to finish our discus- 
sion we would like to show how the obtained polynomials reduce to the usual 
Hermite polynomials in the limit q — > 1~. For this purpose we would like to 
demonstrate that the recursion relation (4.8) is actually the recursion relation 
for the continuous q— Hermite polynomials [42,43]. This means that we have to 
demonstrate that under some conditions (to be specified) the recursion (4.8) is 
equivalent to 

2xH n {x | q) = H n+1 (x | q) + (1 - q n )H n ^(x | q). (4.50) 

known for q-Hcrmitc polynomials. To demonstrate the equivalence we assume 
that x = cos 9 in (4.50) and then, let z — e l9 . Furthermore, we assume that 

H n (x\q)=z n H n (z- 2 ;q), (4.51) 

allowing us to obtain, 

(z + z- x )z n H n = z n+1 H n+1 + (1 - q n )z n - x H n - X . (4.52) 

Finally, we set z _1 = \/t which brings us back to (4.8). This time, however, we 
can use results known in literature for q— Hermite polynomials [41-43] in order 
to obtain at once 

H n (xJ^ | q) 
lim = H n (x), (4.53) 

q^l 1 1-q 

V 2 

where H n (x) are the standard Hermitian polynomials. In view of (4.48), (4.49), 
not surprisingly, the S-W polynomials are also reducible to H n (x). Details can 
be found in the same references. 



4.4 ASEP, q-deformed harmonic oscillator and spin chains 

In this subsection we would like to connect the results obtained thus far with 
the XXX and XXZ spin chains. Although a connection with XXX spin chain 
was established already at the beginning of this section, we would like to arrive 
at the same conclusions using alternative (physically inspired) arguments and 
methods. To understand the logic of our arguments we encourage our readers 
to read Appendix A at this point. In it we provide a self contained summary of 
results related to the asymmetric simple exclusion process (ASEP), especially 
emphasizing its connection with static and dynamic properties of XXX and XXZ 
spin chains. 

ASEP was discussed in high energy physics literature, e.g. see [44], in con- 
nection with random matrix ensembles. To avoid repeats, we would like to use 
the results of Appendix A in order to consider the steady-state regime only. To 
be in accord with literature on ASEP, we would like to complicate matters by 
imposing some nontrivial boundary conditions. 

In the steady -state regime equation (A. 12) of Appendix A acquires the form 
: SC = A. Explicitly, 

SC = p L ED- p R DE. (4.54) 
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In the steady-state regime, the operator S becomes an arbitrary c-number [45] . 
In view of this, following Sasamoto [46] we rewrite (4.54) as 

PrDE — plED = ( (D + E) . (4.55) 

Such operator equation should be supplemented by the boundary conditions 
which are chosen to be as 

a<W\E = {<W\ and [3D \ V >=(\V > . (4.56) 

The normalized steady-state probability for some configuration C can be written 
now as 

p(c) <w\x 1 x 2 ..-x N \v> 
[L) ~ <w\c N \v> [bn 

with the operator Xi being either D or E depending on wether the i— th site is 
occupied or empty. To calculate P(C) we need to determine Q while assuming 
parameters a and (3 to be assigned. We demonstrate in Appendix A that it 
is possible to equate Q to one so that, in agreement with [47], we obtain the 
following representation of D and E operators: 

11 11 

D = - + a, E=- + a+ 

l-q yjl-q l-q yjl-q 

converting equation (4.54) into (4.28d). In view of this mapping into g— deformed 
oscillator algebra, we can expand both vectors | V > and < W | into a Fourier 
series, e.g. | V >= J2 m &m{V) \ m > where, using (4.33), we put | m >= 
By combining equations (4.33) and (4.56) and results of Appendix A we obtain 
the following recurrence equation for fi n : 

Q n (V)(^ - 1) = Q n+1 (V)y/l - q n +\ (4.58) 

Following [47] we assume that < | V >= 1. Then, the above recurrence 
produces 

n n (v) = -_= , 

V(5,9)n 

with parameter v = — 1. Analogously, we obtain: 

n n (w) = 



y/{q,q)n 

with w = — 1. Obtained results exhibit apparently singular behavior for 
q — > 1 _ . These singularities are only apparent since they cancel out when one 
computes quantities of physical interest discussed in both [48] and [49] . As re- 
sults of Appendix A indicate, such a crossover is also nontrivial physically since 
it involves careful treatment of the transition from XXZ to XXX antiferromag- 
netic spin chains. Hence, the results obtained thus far enable us to connect the 
partition function (4.2) (or (4.7)) with either XXX or XXZ spin chains but arc 
not yet sufficient for making an unambiguous choice between these two models. 
This task is accomplished in the rest of this section. 
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4.5 Crossover between the XXZ and XXX spin chains: 
connections with the KPZ and EW equations and the 
lattice Liouville model 



Following Derrida and Malick [49] , we notice that ASEP is the lattice version of 
the famous Kradar-Parisi-Zhang (KPZ) equation [50]. The transition q — ► 1~ 
corresponds to transition (in the sense of renormalization group analysis) from 
the regime of ballistic deposition whose growth is described by the KPZ equa- 
tion to another regime described by the Edwards- Wilkinson (EW) equation. In 
the context of ASEP (that is microscopically) such a transition is discussed in 
detail in [51]. Alternative treatment is given in [49]. The task of obtaining the 
KPZ or EW equations from those describing the ASEP is nontrivial and was 
accomplished only very recently [oliveiraetal, Lazarides]. It is essential for us 
that in doing so the rules of constructing the restricted solid-on -solid (RSOS) 
models were invoked. From the work by Huse [54] it is known that such models 
can be found in four thermodynamic regimes. The crossover from the regime III 
to regime IV is described by the critical exponents of Friedan, Qui and Shcnkcr 
unitary CFT series [55]. The crossover from regime III to regime IV happens 
to be relevant to crossover from the KPZ to EW regime as we would like to 
explain now. 

As results of Appendix A indicate, the truly asymmetric simple exclusion 
process is associated with the XXZ model at the microscopic level and with 
the KPZ equation/model at the macroscopic level. Accordingly, the symmetric 
exclusion process is associated with the XXX model at the microscopic level 
and with the EW equation/model at the macroscopic level. At the level of 
Bethe ansatz for open XXZ chain with boundaries full details of the crossover 
from the KPZ to EW regime were exhaustively worked out only recently [56]. 
For the purposes of this work it is important to notice that for certain values of 
parameters the Hamiltonian of open XXZ spin chain model with boundaries 
can be brought to the following canonical form 



1 1 1 

H XXZ = - [£ (of + ay j+1 + -(q + cT VK+i) + 



•nj 



3 = 1 



(4.59) 

In the case of ASEP we have q — \/pr/pl so that for physical reasons parameter 
q is not complex. However, mathematically, we can allow for q to be complex. In 

7T 

particular, following Pasquer and Saleur [57] we can let q — e 11 with 7 = . 

For such values of q use of finite scaling analysis applied to the spectrum of the 
above defined Hamiltonian produces the central charge 

c=l-_L_ ,^ = 2,3,.... (4.60) 



15 That is equation (1.3) of [56]. 
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of the unitary CFT series. Furthermore, if Ci is the generator of the Tcmpcrlcy- 
Lieb algebra 16 !, then Hxxz can be rewritten as [58] 

N-l 

H XX z = -J2[e i --(q + q- 1 )}. ( 461 ) 
i=i 

This fact allows us to make immediate connections with quantum groups and 
theory of knots and links. Below, in Section 5 we shall use different arguments to 
arrive at similar conclusions. The results just described allow us to connect the 
CFT and exactly integrable lattice models. If this is the case, one can pose the 
following question: given the connection we just described, can we write down 
explicitly the corresponding path integral string-theoretic models reproducing 
results of exactly integrable lattice models at and away from criticality? Before 
providing the answer in the following subsection, we would like to conclude 
this subsection with a partial answer. In particular, we would like to mention 
the work by Faddeev and Tirkkonen [59] connecting the lattice Liouville model 
with the spin 1/2 XXZ chain. Based on this result, it should be clear that in 
the region c < 1 it is indeed possible by using combinatorial analysis described 
above to make a link between the continuum and the discrete Liouville theo- 
rie^Zl- It can be made in such a way that, at least at crtiticality, the results 
of exactly integrable 2 dimensional models are in agreement with those which 
are obtainable field- theoretically. The domain c > 1 is physically meaningless 
because the models (other than string-theoretic) we discussed in this section 
loose their physical meaning in this region. This conclusion will be further 
reinforced in the next subsection. 

4.6 ASEP, vicious random walkers and string models 

We have discussed at length the role of vicious random walkers in derivation 
of the Kontsevich-Witten (K-W) model in our previous work [60]. Forrester 
[61] noticed that the random turns vicious walkers model is just a special case 
of ASEP. Further details on connections between the ASEP, vicious walkers, 
KPZ and random matrix theory can be found in the paper by Sasamoto [62]. 
In the paper by Mukhi [63] it is emphasized that while the K-W model is the 
matrix model representing c < 1 bosonic string, the Penner matrix model 
with imaginary coupling constant is representing c = 1 Euclidean string on 
the cylinder of (self-dual) radius R = l|^f| Furthermore, Ghoshal and Vafa [65] 
have demonstrated that c — 1 , R = 1 string is dual to the topological string on 
a conifold singularity. We shall briefly discuss this connection below. Before 
doing so, it is instructive to discuss the crossover from c = 1 to c < 1 string 
models in terms of vicious walkers. To do so we shall use some results from our 
work on K-W model and from the paper by Forrester [61]. 

ie That is e| = e^, e^e^iei = q~ 1 ei and eiej = eja for \i — j\ > 2. 
17 The matrix c = 1 theories will be discussed separately below. 
18 This was initially demonstrated by Distler and Vafa [64]. 
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Thus, we would like to consider planar lattice where at the beginning we 
place only one directed path P: from (a, 1) to (b, Ny^ The information about 
this path can be encoded into multiset Hor y (P) of y-coordinates of the hori- 
zontal steps of P. Let now 

w(P)= 11 x t . (4.62) 

i=Hor y (P) 

Using these definitions, the extension of these results to an assembly of directed 
random vicious walkers is given as a product: W(P) = w(Pi) ■ ■ ■ w(Pk). Finally, 
the generating function for an assembly of such walkers is given by 

h b - a ( Xl ,...,x N )=J2 W (P), ( 4 - 63 ) 
P 

where W(P) is made of monomials of the type N provided that 

mi + • • • + rrcjv = b — a. The following theorem [ , ] is of central importance for 
calculation of such defined generating function. 

Given integers < a\ < ■ ■ ■ < au and < b\ < ■ ■ ■ < bk , let M^j be the 
k x k matrix Mjj = -04(^11 ■■■i x n) then, 

detM = W(P) (4.64) 
P 

where the sum is taken over all sequences (Pi, ...,Pfc) = P of nonintersecting 
lattice paths Pi : (<ij,l) — * (bi,N),i = 1 — k. 

Let now a% — i and bj — Aj + j so that 1 < i,j • < k with A being a 
partition of N with k parts then, detM = s\(xi, Xn), where s\(x) is the 
Schur polynomial. In our work [60 ] we demonstrated that in the limit N — > 00 
such defined Schur polynomial coincides with the partition(generating) function 
for the Kontsevich model. Many additional useful results related to Schur 
functions are discussed in our recent paper [2]. 

To get results by Forrester requires us to apply some additional efforts. 
These are worth discussing. Unlike the K-W case, this time, we need to discuss 
the continuous random walks in the plane. Let x-coordinate represent "space" 
while y-coordinate- "time". If initially (t = 0) we had k- walkers in the positions 
— L < x\ < X2 < ■ ■ ■ < Xk < L, the same order should persist V t > 0. 
At each tick of the clock each walker is moving either to the right or to the 
left with equal probability p (that is we are in the regime appropriate for the 
XXX spin chain in the ASEP terminology). As before, let xo = (xi t o, ...,Xk,o) 
be the initial configuration of A:— walkers and xy = [xxj, ...,Xkj) be the final 
configuration at time t. To calculate the total number of walks starting at t = 
at x and ending at time t at Xj we need to know the probability distribution 
Wk(x.Q — + x /ji) that the walkers proceed without bumping into each other. 

19 Very much in the same way as discussed already in Section 2. 
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Should these random walks be totally uncorrelated, we would obtain for the 
probability distribution the standard Gaussian result: 



Wk{X ° X/ '' } = ' (465) 

In the present case the walks are restricted (correlated) so that the probability 
should be modified. This modification can be found in the work by Fisher and 
Huse [66] . These authors obtain 



W k (x ^Xf,t) = U k {x ,x f ;t) J^J," ( 466 ) 



x exp{-( X 2 +xl)/2Dt} 
with 

tf fc (xo,x /; t) = E £(g)exp[ (x/ ^ Xo) ]. (4.67) 

ges k 

In this expression s(g) = ±1, and the index g runs over all members of the 
symmetric group S k . Mathematically, following Gaudin [67], this problem can 
be looked upon as a problem of a random walk inside the k— dimensional kalei- 
doscope (Weyl cone) usually complicated by imposition of some boundary con- 
ditions at the walls of the cone. Connection of such random walk problem with 
random matrices was discussed by Grabiner [68] whose results were very re- 
cently improved and generalized by Krattenthallcr [69] . In the work by de Haro 
some applications of Grabiner's results to high energy physics were considered 
[70]. Here we would like to approach the same class of problems based on the 
results obtained in this paper. In particular, some calculations made in [66] 
indicate that for L — > oo with accuracy up to O (L 2 / Dt) it is possible to rewrite 
f/fc(xo,x/;£) as follows: 

f7fe(xo,x/;i) ~ consfA( X/ )A(x )/ (£>£)"* +0(L 2 /Dt) (4.68) 

with n k = (l/2)fc(fc - 1) and const = 1/112! ■ • • (k - 1)! and A(x) being the 
Vandermonde determinant, i.e. 

A(x) = Y[( Xi - Xj ). (4.69) 

i<j 

Next, from standard texts in probability theory it is known that non-normalized 
expression, say for W k (xo — » x/ ; t) in the limit of long times provides the 
number of random walks of n steps (since n t) from point x to point x/. 
Hence, the same must be true for Wfc(xo — ► x/;t) and, therefore, Wfe(xo — > 
x/; £) « det M. Consider such walks for which x(i = 0) = xo = x(t = tf) = xy. 
Then, using (4.66) and (4.68) we obtain the probability distribution for the 
Gaussian unitary ensemble [71], i.e. 

W fc (x = xf,t) = const A 2 (x) exp(-x 2 ). (4.70) 
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Some additional manipulations (described in our work [60] ) using this ensemble 
lead directly to the K-W matrix model. Forrester [61] had considered a related 
quantity: the probability that all k vicious walkers will survive at time tf. To 
obtain this probability requires integration of Wfefxp — > x.f,t) over the simplex 
A defined by — L < x\ < x-i < ■ ■ ■ < x^ < lF°l Without loss of generality, it 
is permissible to use Wt(x — x/;t) instead of Wk(*o — ► x /!^) m calculating 
such a probability. Then, the obtained result coincides (up to a constant) with 
the partition function of topological gravity Z(g), equation (3.1) of [72,73^]]). 
Furthermore, such defined partition function can be employed to reproduce back 
the Hermite polynomial Hk(x) defined by (4.53) which has an interpretation as 
the wavefunction(amplitude) of the FZZT D— brane [72,73]. Specifically, we 
have 

< det(x - M) >= (|) 2 H k {x^) 

= tt-I-t fdM detfx- M)e-3 trM \ (4.71) 
z \9) J 

This expression is a special case of Heine's formula representing monic orthog- 
onal polynomials through random matrices. In the above formula k is related 
to the size of Hcrmitian matrix M and g is the coupling constant. 

Following Forrester [61], the result (4.66) can be treated more accurately 
(albeit a bit speculatively) if, in addition to the parameter D we introduce 
another parameter a - the spacing between random walkers at time t = 0. 
Furthermore, if the time direction is treated as space direction (as it is commonly 
done for Id quantum systems in connection with 2d classical systems), then yet 
another parameter r(k,t) should be introduced which effectively renormalizcs 
D. This eventually causes us to replace Z(g) by the following integral (up to a 
constant) 

k 00 

Z{g) = ft / <bi exp(-^ In 2 x t ) Y[( Xi - x 3 f = f dMe~^ tr{ln M)2 (4.72) 

Tierz demonstrated [74] that Z(g) (up to a constant) is partition function of the 
Chern-Simons (C-S) field theory with gauge group U(k) living on the 3-sphere 
S 3 . Okuyama [73] used (4.71) in order to get analogous result for a D-brane am- 
plitude in C-S model. Using Heine's formula, he obtained the Stieltjes-Wiegert 
(S-W) polynomial, our equation (4.47), which can be expressed via the Rogers- 
Szego polynomial according to (4.49) and, hence, via the q— Hermite polyno- 
mial in view of the relation (4.51). Since in the limit q — > l _ the q— Hermite 



20 Such type of integration is described in detail in our papers, Parts I and II, from which 
it follows that in the limit L — > oo such a simplex integration can be replaced by the usual 

fc oo fc 

integration, i.e Y[^ x i ' " ^ TT / n^ x '» ' " ' n accord with Forrester. 

i —aoi 

21 Since the hermitian matrix model given by (3.1) is just a partition function for the Gaus- 
sian unitary ensemble [71]. 
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polynomial is reducible to the usual Hermite polynomial according to (4.53), 
there should be analogous procedure in going from the partition function Z(g) 
to Z{g). Such a procedure can be developed, in principle, by reversing argu- 
ments of Forrester. However, these arguments are much less rigorous and phys- 
ically transparent than those used in previous subsection where we discussed 
the crossover from XXZ to XXX model. In view of the results presented in the 
following section, we leave the problem of crossover between the matrix ensem- 
bles outside the scope of this paper. To avoid duplications, we refer our readers 
to the paper by Okuyama [73] where details are provided relating our results to 
the topological A and B -branes. 

5 Gaudin model as linkage between the WZNW 
model and K-Z equations. Recovery of the 
Veneziano-like amplitudes 



5.1 General remarks 

We would like to remind to our readers that all results obtained thus far can be 
traced back to our equation (4.7) defining the Rogers-Szego polynomial which 
physically was interpreted as partition function for the ferromagnetic P-F spin 
chairF^I. In previous sections numerous attempts were made to connect this 
partition function to various models, even though already in Section 4.1 we came 
to the conclusion that in the limit of infinitely long chains the antiferromagnetic 
version of P-F spin chain can be replaced by the spin 1/2 antiferromagnetic 
XXX chain. If this is so, then fom literature it is known that behaviour of such 
spin chain is described by the SU\{2) WZNW model [30]. Hence, at the physical 
level of rigor the problem of connecting Veneziano amplitudes to physical model 
can be considered as solved. In this section we argue that at the mathematical 
level of rigor this is not quite so yet. This conclusion concerns not only problems 
dicussed in this paper but, in general, the connection between the WZNW 
models, spin chains and K-Z equations. It is true that K-Z equations and 
WZNW model are inseparable from each other [30] but the extent to which spin 
chains can be directly linked to both the WZNW models and K-Z equations 
remains to be investigated. We would like to do so in this section. For the 
sake of space, we shall discuss only the most essential facts leaving (with few 
exceptions) many details and proofs to literature. 

Following Varchenko [8] , we notice that the link between the K-Z equations 
and WZNW models can be made only with help of the Gaudin model, while 

22 The antiferromagnetic version of P-F spin chain is easily obtainable from this ferromag- 
netic version as discussed in Section 3. 
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the connection with spin chains can be made only by using the quantum ver- 
sion of the K-Z equation. Such quantized version of the K-Z equation is not 
immediately connected with the standard WZNW model as discussed in many 
places [8,75]. In this section, we would like to discuss in some detail the Gaudin 
model and its relation to the P-F spin chain and, hence to the Veneziano model 
formulated in Part II. We begin with summary of facts related to this model. 

5.2 Gaudin magnets, K-Z equation and P-F spin chain 

Although theory of the Gaudin magnets plays an important role in topics such 
as Langlands correspondence, Hitchin systems, etc. [76-78] in this work we do 
not discuss these topics. Instead, we would like to focus only on issues of 
immediate relevance to this paper. Gaudin came up with his magnetic chain 
model in 1976 [67] being influenced by earlier works of Richardson [79, 80] on 
exact solution of the BCS equations of superconductivity. This connection with 
superconductivity will play an important role in what follows. 

In physics literature all Gaudin-type models are based on the SU(2) algebra 
of spin operatorj^l. Instead of one Hamiltonian, the set of commuting Hamil- 
tonians of the type [82] 

N 3 
j&i)=l <x=l 

is used. In view of the fact that, by construction, [Hi,Hj] = 0, 3N(N — 1), the 
coefficients wfj should satisfy the following equations 

<j w ]k + w % w lk ~ w ik™% = 0. (5.2) 

These equations can be solved by imposing the antisymmetry requirement: 
w°j = — which can be satisfied by replacing wfj by the unknown func- 
tions wfj — f a (zi — Zj) of difference between two new real parameters Zi 
and Zj. It is only natural to make further restrictions based on requirement 
that the z— component of the total spin S 3 =XL °"f is conserved. This causes 
wjj = wfj = and wfj = Yij thus leading to equations 

Y^Xjk + Y^Xjk + X^Xij = 0. (5-3) 

These constraint equations admit the following sets of solutions: 

X^ = Y^ = (rational), (5.4a) 

Zi Zj 

Xij = — — , Y^ = cos (zi — Zj) (trigonometric), (5.4b) 



23 In mathematics literature to be used below [8,75] the SL(2, C) group is used instead of 
its subgroup, SU(2) [ 81]. 
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Xij = r-ry- r> Yij = cosh (2^ — z^) (hyperbolic). (5.4c) 

Slllll yZ>i j ) 

While the first solution, (5.4a), to be used in this work, corresponds to the long 
range analog of the standard XXX spin chain, the remaining two solutions 
correspond to the long range analog of the XX Z spin chain. 

Folloving Varchenko [8] we are now in the position to write down the K-Z 
equation. For this purpose we combine equations (5.1) and (5.4a) and reintro- 
duce the coupling constant g (so that wfj — > gw"A in such a way that the K-Z 
equation acquires the form 

d 

(k- Hifa,..., z N ))Hzi,..., zn) = 0, i = l,...,N,. (5.5) 

ozi 

where k = g . This result requires several comments. First, from the theory 
of WZNW models it is known that parameter k cannot take arbitrary values. 
For instance, for SZ7i(2) WZNW model k = § [30]. Second, we can always 
rescale z-coordinates and to redefine the Hamiltonian to make the constant 
arbitrary small. Apparently, this was asssumed in the asymptotic analysis of 
the K-Z equation described in [7,8]. Third, if this is the case, then such analysis 
(to be used below) differs essentially from other approaches connecting string 
models with spin chains, e.g. see [83], because such a connection was made in 
these works for SU (N)-type magnets (or gauge theories) in the unphysical limit 
N — > oo. Since for SU(N) models n = + N), in the limit N ~ + oo we have 
k — ► oo . The WKB-type method of Rcshctikhin and Varchenko (to be discussed 
below) fails exactly in this limit. 

With K-Z equation defined, we would like to make a connection between 
the Gaudin and P — F model. To a large extent this was already accomplished 
in [84]. Following this reference, we define the spin Calogero (S-C) model as 
follows 

H s-c = ~ Ytf + ^ + 9 E 7^Z% (5-6) 

i i<j v-i z 3) 

to be compared with TL in (3.20. Using the rational form of the Gaudin 
Hamiltonian, this result can be equivalently rewritten as 

Hs-c = \ + ^xf + i^lpuHt}). (5.7) 

i 

That this is indeed the case can be seen by the following chain of arguments. 

Consider the strong coupling limit (g — > oo) of Hs-c so that the kinetic 
term is a perturbation. Next, we consider the eigenvalue problem for one of the 
Gaudin's Hamiltonians, i.e. 

Hl y(l) = E®i&W, (5.8) 



24 We added the oscillator-type potential absent in the original work [84] for the sake of 
additional comparisons, e.g. with (3.4). In what follows such a constraint is not essential and 
will be ignored. 
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and apply the operator ipi to both sides of this equation. Furthermore, consider 
in this limit the combination Hs~c^^- Provided that the eigenvalue problem 
(5.8) does have a solution, it is always possible to Fourier expand (ipi^^) using 
as basis set In such a case we end up with the eigenvalue problem for the 
P-F spin chain in which the eigcnfunctions are the same as for the Gaudin's 
problem and the eigenvalues are ipiE^ . Physical significance of this result will 
be discussed in detail below. Before doing so, we have to make a connection 
between the K-Z equation (5.5) and Gaudin eigenvalue equation (5.8) following 
[7,8]. 

We begin by replacing SU(2) spin operators by SL(2, C) = sl 2 operators e, 
/ and h obeying following commutation relations 

[h,e]=2e; [e,f]=h; [h,f] = -2f. (5.9) 

This Lie algebra was discussed in our previous work, Part II, in connection with 
new models reproducing Veneziano amplitudes. In this work, we shall extend 
these results following ideas of Richardson and Varchenko. 

From [81] it is known that SU(2) is just a subgroup of sl 2 - Introduce the 
Casimir element Q £ sl 2 ® sl 2 via 

= e®/ + /Oe+^/i(g)/i (5.10) 

so that Vx € sl 2 it satisfies the commutation relation [x £g> 1 + 1 (g> x, f2] =0 
inside theU(sl2) <S> U{sl 2 ) where U(sl 2 ) 1S the universal enveloping algebra of 
sl 2 . Consider the vector space V = V\ ® V 2 ® • • • <8> Vat. An element x <G sl 2 
acts on V as follows: a;®l®l(g)---®l + -- - + l®l®---®a;. For indices 
1 < i < j < N let fiW) : V — > V be an operator which acts as f2 on i-th and 
j-th position and as identity on all others, then the K-Z equation can be written 
as 

« * = £ , = 1,...,^. (5.11) 

In the simplest case, the K-Z equation is defined in the domain U — {(z\, zat) e 
C w I Zj }. 

From now on we shall use equation (5.11) instead of (5.5). To connect K-Z 
equation with the XXX Gaudin magnet we shall use a kind of WKB method 
developed by Reshetikhin and Varchenko [7] and summarized in lecture notes 
by Varchenko [8]. Following these authors, we shall look for a solution of (5.11) 
in the form (k —> 0) : 

$(z,«) = e^ z >{/o(z) + k/i(z) + •••}, (5.12) 

where z = {z\, Zn}, S(z) is some scalar function (to be described below) 
and fj(z), j = 0, 1,2, are V— valued functions. Provided that the function 
S is known, V— valued functions can be recursively determined (as it is done in 



30 



WKB analysis). Specifically, given that Hi — j we obtain, 

jjti z i — z j 

ffi/oW = Q^foW (5-13) 
to be compared with (5.8). Next we get 

and so on. Since the function S(z) (the Shapovalov form) plays an important 
role in these calculations, we would like to discuss it in some detail now. 

5.3 The Shapovalov form 

Let us consider the following auxiliary problem. Let A(x) and B(x) be some 
pre assigned polynomials of degree n and n — 1 respectively. Find a polynomial 
C(x) of degree n — 2 such that the differential equation 

A(x)y" (x) - B{x)y'{x) + C{x)y{x) = (5.15) 

has solution which is polynomial of preassigned degree k. Such polynomial 
solution is called the Lame' function. Stieltjes [7,8] proved the following 

Theorem. Let A and B be given polynomials of degree n and n — 1, re- 

fn ■ 

spectively so that B(x)/A(x) = Y^j=i ~- Then there is a polynomial C of 

degree n — 2 and a polynomial solution y(x) — Y\^ =1 (x — Xi) of (5.15) if and 
only if x =(x\, xu) is the critical point of the function 

k n 

<S> Kn (x 1 ,...,x k ;z 1 ,...,z n ) = \\\\{x J -z l )- m > \\ (x.-xtf. (5.16) 

j=l »=1 l<i<j<k 

A point x is critical for $(x) if all its first derivatives vanish at it. 

We would like now to make a connection between the Shapovalov form S and 
the results just obtained. S is symmetric bilinear form on previously introduced 
space V such that S(v,v) — 1, S(hx,y) = S(x,hy), S(ex,y) — S(x,fy) where 

h, e,f are defined in (5.9). Furthermore, S(fl(x\ ® X2), y\ ® 2/2) = S{x\ ® 
X2, ^(yi (E> 2/2)) Vzi, yi e Vi and Vx 2 , y2 € V2. As result, we obtain, 

S(H i x,y) = S(x,H i y)\fx,yeV. (5.17) 

Next, let m be some nonnegative integer and V m be the irreducible Verma 
module with the highest weight m and the highest weight singular vector v m , 

i. e. 

hv m = mv m , ev m = 0. (5.18) 
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Consider a tensor product V = V® M — V mi ®---®V mn so that M — (mi, m n ). 
VV mi vectors v mi , fv mi , f 2 v mi ,- ■ ■,f nii v mi form a basis of V mi r 5 bo that the 
Shapovalov form is orthogonal with respect to such a basis and is decomposable 
as S = S mi <8> • • • <8> S mn Let, furthermore, J = (ji,—,j n ) be a set of nonnegative 
integers such that ji + ■ ■ ■ + j n = k where k is the same as in (5.16) and 
< ji < mi- This allows us to define the vectors / j vm — f^v mi ® • • • ® P n v mn 
These vectors {J j Vm} are by construction orthogonal with respect to the 
Shapovalov form and provide a basis for the space V® M . Introduce the weight 
of a partition A as \A\ = a\ + a 2 + ... then, in view of (5.18), we define the 
singular vector f J VM via 

h(f J v M ) = (\M\ - 2 \J\)f J v M , e(fv M ) = (5.19) 

of weight \M\ — 2 | J| F^l The Bethe ansatz vectors V for the Gaudin model can 
be defined now as 

V(x,z) = ^A/(x,z)/ J « M (5.20) 
,/ 

where x is a critical point of $(x, z) which was defined by (5.16). A function 
Aj(x, t) is defined as follows 

k 

Aj(Z,z)= tl x .-\ . ( 5 - 21 ) 

aeV(k;J) 4=1 Xl Z<T W 

with Vik\ J) being the set of maps a from the {1, k} to {1, n}. Finally, 
using these definitions it is possible to prove that 

d 2 

S(V(x,z), V(x,z)) = det ( ln$ ktn (xi,...,x k ;zi,...,z n )). (5.22) 

l<ij<fe OXiOXj 

The equations determining critical points 

^^0)^:^(^)^)1^0=0 (5.23) 

are the Bethe ansatz equations for the Gaudin model. Using these equations 
the eigenvalue equation (5.13) for the Gaudin model now acquires the following 
form 

i/ J (z°)V(x°,z ) = JLln<Z> k , n (x(z), z) | z=z o V(x°,z°). (5.24) 

In the next subsection we shall sudy in some detail the Bethe ansatz equation 
(5. 23). This will allow us to define eigenvalues in (5.24) explicitly. 



25 According to [8] in all subsequent calculations it is sufficient to use the finite Verma 
module, i.e. L m = V m / < f m+1 Vm > ■ This restriction is in accord with our previous 
calculations, e.g. see Part II, Section 8, where such a restriction originates from the Lefschetz 
isomorphism theorem used in conjunction with supersymmetric model reproducing Veneziano 
amplitudes. 

2e This fact can be easily undestood from the properties of sl2 Lie algebra representations 
since it is known, [8] and Part II, that for the module of highest weight m we have h(f k v m ) = 
(m - 2k)(f k v m ). 
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5.4 Mathematics and physics of the Bethe ansatz equa- 
tions for XXX Gaudin model according to works by 
Richardson. Connections with the Veneziano model 

Using (5.16) in (5.23) produces the following set of the Bethe ansatz equations: 

n k „ 

£^^=E^-' j = l,..;k. (5.25) 

i=l 3 1 i=l 3 1 

To understand the physical meaning of these equations we shall use extensively 
results of two key papers by Richardson [79,80]. To avoid duplications, and for 
the sake of space, our readers are encoraged to read thoroughly these papers. 
Although originally they were written having applications to nuclear physics in 
mind, they are no less significant for condensed matter [82] and atomic physcs 
[85] . Because of this, only nuclear physics terminology will be occasionally used. 
At the time of writing of these papers, QCD was still in its infancy. Accord- 
ingly, no attempts were made to apply Richardson's results to QCD. Recently, 
Ovchinnikov [86] have conjectured that the Richardson-gaudin equations can be 
useful for development of color superconductivity in QCD [87] . Incidentally, in 
the same paper [87] it is emphasized that such type of superconductivity can 
exist only if the number of colors is not too large, e.g. N c =3. This fact is in ac- 
cord with remarks made in Section 5.2 regarding the validity of the WKB-type 
methods in the limit N— > oo for the K-Z equation. 

Thus, following Richardson [80], we consider the system of interacting 
bosons described by the (pairing) Hamiltoniarf^l 

Here we have hi — a k ak ' = S a k a ^k an< ^ — Yl a -kflk- 

k(e k =ej) k(e k =ei) k(e k =£() 

It is assumed that the single-particle spectrum {ei} is such that Si < £;+i V7 
and that the degeneracy of Z-tli level is fi; so that the sums (over k ) each 
contain f2j terms. It is assumed furthermore that the system possesses the 
time-reversal symmetry implying £k = £-k- The operators and <Zk obey 
usual commutation rules for bosons, i.e. [ak,»j^] = ^kk'- The sign of the 
coupling constant in principle can be both positive and negative. We shall 
work, however, with more physically interesting case of negative coupling (so 
that g in (5.26) is actually \g\). 

An easy computation using commutation rule for bosons produces the fol- 
lowing results 

[h l ,A+}=2S ll ,A+, (5.27a) 
[A h A+]=2d w (Q l +2h l ), (5.27b) 

27 In the paper with Sherman [79] Richardson explains in detail how one can map the 
fcrmionic (pairing) system into bosonic. 
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[n l ,A l ,]=-26 w A l . 



(5.27c) 



If we make a replacement of hi in (5.27a) and (5.27.c) by ^ + hi = — and 

keep the same notation in the r.h.s. of (5.27b) we shall arrive at the SI2 Lie 
algebra isomorphic to that given in (5.9). The same Lie algebra was uncovered 
and used in our Part II for description of new models describing Veneziano 
amplitudes. Because of this, we would like now to demonstrate that the rest 
of arguments of Part II can be implemented now in the present context thus 
making the P-F model (which is derivative of the Richardson-Gaudin XXX 
model) correct model related to Veneziano amplitudes. 

Following Richardson [80], we notice that the model described by Hamil- 
tonian (5.26) and algebra (5.27) admit two types of excitations: those which 
are associated with the unpaired particles and those with coupled pairs. The 
unpaired ^—particle state is defined by the following two equations 

h\<p v >=v\<p v >, (5.28) 

Ai I <p v >= V/. (5.29) 
Here, h = J2i ™i so that, in fact, 

hi I ip v >= vi\Lp u > (5.30) 

and, therefore, v = v\. Furthermore 

H\y v >=Y Jl zin\v„> ■ (5-31) 

Following Richardson, we want to demonstrate that parameters Ei in (5.31) can 
be identified with parameters zi in the Bethe equations (5.25). Because of this, 
the eigenvalues for the P-F chain are obtained as described in Section 5.2., that 
is 

These are eigenvalues of hi defined in (5.30). Furthermore, this eigenvalue 
equation is exactly the same as was used in Part II, Section 8, with purpose 
of reproducing Veneziano amplitudes. Moreover, equations (5.28) and (5.29) 
have the same mathematical meaning as equations (5.19) defining the Verma 
module. Because of this, we follow Richardson's paper to describe this module 
in physical terms. By doing so additional comparisons will be made between 
the results of Part II and works by Richardson. Since the Hamiltonian (5.26) 
describes two kinds of particles: a) pairs of particles (whose total linear and 
angular momentum is zero) and, b) unpaired particles (that is single particles 
which do not interact with just described pairs), the total number of (quasi) 
particles is n = N + Since we redefined the number operator as ^ + hi = 

28 In Richardson's paper we find instead: n = 2N + v. This is, most likely, a misprint as 
explained in the text. 
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= N; we expect that , once the correct state vector describing excitations 

is found, equation (5.30) should be replaced by the analogous equation for N; 
whose eigenvalues will be % + ^>F^I 

A simple minded way of creating such a state is by constructing the following 
state vector ■ ■ ■ Af N \ ip v > . This vector does not possess the needed 
symmetry of the problem. To create the state vector (actually, the Bcthc vector 
of the type given by (5.20)) of correct symmetry one should introduce a linear 
combination of At operators according to the following prescription: 



K=Y. ««(0^i"> « = 1. -> N ( 5 - 33 ) 

i 



with constants u a (l) to be determined below. The (unnormalized) Bethe-type 
vectors are given then as | ip >= ■ ■ ■ B^ \ <p v > and, accordingly, instead of 
(5.31), we obtain 

H I V> >= E, I ^ > + l H > B t ■ ■ ■ B n] I Vu > • (5-34) 

The task now lies in calculating the commutator and to determine the constants 
u a (l). Details can be found in Richardson's paper [80]. The final result looks 
as follows 

H | tp>-E\ip> (5.35) 

N 

= E(Il i? 7)EA + [( 2 £i- £; «K(o+E( r!i ' + 2 ^') u «( r )+4g Yl Mp a \\v v >. 

a=l 7 ^q I I' P(P^a) 

By requiring the r.h.s. of this equation to be zero we arrive at the eigenvalue 
equation 

JV 

H | V >= E | V >, where E = + ^ E a- (5.36) 

l a = l 

Furthermore, this requirement after several manipulations leads us to the Bethe 
ansatz equation^] 

2~g + E E^E~ a -^2^E: = ' a=1 — N > (5 - 37a) 

y /3(/3jbat) 1 1=1 

as well to the explicit form of coefficients u a (l) : u a (l) = l/(2e; — E a ) and 
the matrix elements M a _p (since, by construction, u a (l)up(l) = M a> fjU a {l) + 



29 These amendments are not present in Richardson's paper but they are in accord with its 
content. 

30 It should be noted that in the original paper [80] the sign in front of the 3rd term in 
the l.h.s. is positive. This is because Richardson treats both positive and negative couplings 
simultaneously. Equation (5.37a) is in agreement with (3.24) of Richardson-Sherman paper 
[79] where the case of negative coupling (pairing) is treated. 
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Mp j0l up(l)). In the limit g — > we expect E a — > 2e; and — > in accord with 
(5.28)-(5.30). Therefore, we conclude that %-+z^ is an eigenvalue of the operator 
N; acting on | tp > in accord with remarks made before. In the opposite limit: 
g — ► oo the system of equations (5.37a) will coincide with (5.25) upon obvious 
identifications: x a E ai 2si +± zi,N k,L +± n and Slz/2 + ^ m;. Next, 
in view of (5.32) and (5.36) we obtain the following result for the occupation 
numbers: 

d 



9£j . 

Z a— 1 



7 ^ 
„ i ^ ' 



(5.38) 



Based on the results just obtained, it should be clear that, actually, e\ ' 

JV 

=i/j + ^ so that ^ = ^ ^ 2L - Richardson [jmp] cleverly demonstrated that 

N 

the combination 7^ must be an integer. 

a=l 

Consider now a special case: N = 1. Evidently, for this case, the derivative 
should also be an integer. For different e-s these may, in general, be 
different integers. This fact has some physical significance to be explained below. 

To simplify matters, by analogy with theory of superconducting grains [82] , 
we assume that the energy £j can be written as £j = d(2i — L— 1), i = 1, 2, i. 
The adjustable parameter <i measures the level spacing for the unpaired par- 
ticles in the limit g — > 0. With such simplification, we obtain the following 
BCS-type equation using (5.37) (for N = 1): 



where G is the rescaled coupling constant. Such an equation was discussed 
in the seminal paper by Cooper [88] which paved a way to the BCS theory of 

L 

superconductivity. To solve this equation, let now F(E) — ^ ^i(2ei — E)^ 1 so 

i=i 

that (5.39) is reduced to 

F(E)=G-\ (5.40) 

This equation can be solved graphically as depicted below 

As can be seen from Fig.l, solutions to this equation for G = oo can be 
read off from the x axis. In addition, if needed, for any N > 1 the system of 
equations (5.37a) can be rewritten in a similar BCS-like form if we introduce 
the renormalized coupling constant G a via 

N 1 

G a = G{l + 2G £ ^— (5.41) 
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Figure 1: Graphical solution of the equation (5.40) 



so that now we obtain: 

F(E a )=Gz\a=l,...,N. (5.37b) 

This sustem of equations can be solved iteratively, beginning with equation 
(5.40). There is, however, better way of obtaing these solutions. In view of 
equations (5.15), (5.16) and (5.23) solutions {E a } of (5.37.b) are the roots of 
the Lame'— type function which is obtained as solution of (5.15). Surprisingly, 
this fact known to mathematicians for a long time has been recognized in nuclear 
physics literature only very recently [89] . 

5.5 Emergence of the Veneziano-like amplitudes as con- 
sistency condition for N — 1 solutions of the K-Z equa- 
tions. Recovery of the pion-pion scattering amplitude 

Since results for the Richardson-Gaudin (R-G) model are obtainable from the 
corresponding solutions of the K-Z equations in this subsection we would like 
to explain why ./V = 1 solution of the Bethe-Richardon equations can be linked 
with the Veneziano-like amplitudes describing the pion-pion scattering. In doing 
so, we shall by pass the P-F model since, anyway, it is obtainable from the R-G 
model. 

Thus, we begin again with equations (5.10)-(5.11). We would like to look 
at the special class of solutions of (5.11) for which the parameter | J| in Verma 
module (5.19) is equal to one. This corresponds exactly to the case N = 1. 
Folloving Varchenko [8], by analogy with (5.16) we introduce the function $(z, t) 
via 

$(z,t)= 11 (zi-Zj) k Y[(t-zi)~K. (5.42) 

\<i<j<L 1=1 
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It is a multivalued function at points of its singularities at z\,...,Zl- Using this 
function, we define the set of 1-forms via 

Wj = $(z,t)-^-, j = l,...,L, (5.43) 

t - Zj 

and the vector of integrals 1^ = (Ii, = (J 7 wi, / Wi) with 7 being 

a particular Pochhammer countour: a double loop winding around any two 
points z a , zp taken from the set Zi, Zl- Deatails can be found in [8,75]. 

We want now to design the singular Verma module for the K-Z equations 
using equation (5.19) and results just presented. Taking into account the 
following known relations: 

a) ef k v m = k(m-k + l)/* -1 ^, and b) hfv m = (m - 2k)f k v m 

for the Lie algebra sh also used in Part II, Section 8, and taking into account 



that in the present (N — 1) case the basis vectors f' 1 'vm — P x v mi 
acquire the following form: f 1 VM = v mi (g> ■ ■ ■ (g> fv ms <g> ■ • ■ ® v mn , s = 1, L, 
provided that m^s are the same as in (5.25) (or (5.42)), the singular vector for 
such a Verma module is given by 

L 

w(l) = IsVm i ® " " " ® K, ( 5 - 44 ) 

8=1 

In view of the Lie algebra relations just introduced, we obtain e • w = or, 
explicitly, 

L 

J2 m sls^0. (5.45) 

s=l 

Hence, for a fixed Pochhammer contour 7 there are L — 1 independent basis 
vectors {w 2 }. They represent L — 1 independent solutions of the K-Z equation 
of the type k = 1 (or iV = 1). Let now z^s be ordered in such a way that z\ < 
■■■ < zl. Furthermore, in view their physical interpretation described in previous 
section, these z[s can be chosen to be equidistant. Consider then a special set 
of Pochhamer contours {7,} around points Zi and i = 1,2, L — 1 and 
consider the matrix M made of integrals of the type MJ = j Uj then, 

any (k = 1)— type solution (f) l (i = 1,2, ...,L — 1) of the K-Z equation can be 
represented as 

^ =^2M l j w j , i = l,2,...,L-l. (5.46) 

3 

From linear algebra it is known that in order for these K-Z solutions to be 
independent we have to require that det M/0. The proof of this fact is given 
in Appendix B. Calculation of the determinant of M is described in detail in [8] 
with the result: 

det M = ±A-i ^ r^rrj s_i 5.47 

r(i-JM) 
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with ±A being some known constantPl and T(x) being Euler's gamma function. 
For L = 2 without loss of generality one can choose z\ = and Z2 = 1, then in 
thus obtained determinant one easily can recognize the Veneziano-type 7r + 7r~ 
scattering amplitude used in the work by Lovelace [90] . We have discussed this 
amplitude previously in connection with mirror symmetry issues [91]. This time, 
however, we would like to discuss other topics. 

In particular, we notice first that all mesons are made of two quarks. Specif- 
ically, we have ud for 7r + , du for ir~ and dd for ir°. These are very much like the 
Cooper pairs with qq quark pairs contributing to the Bose condensate which 
was created as result of spontaneous chiral symmetry breaking. As in the case 
of more familiar Bose condensate, in addition to the ground state we expect to 
have a tower of the excited states made of such quark pairs. Experimentally, 
these are interpreted as more massive mesons. Such excitations are ordered 
by their energies, angular momentum and, perhaps, by other quantum num- 
bers which can be taken into account if needed. Color confinement postulate 
makes such a tower infinite. Evidently, the Richardson- Gaudin (R-G) model fits 
ideally this qualitative picture. Equation (5.40) describes excitations of such 
Cooper-like pairs (even in the limit: G — > 00) as can be seen from Fig.l. In the 
P-F model the factor f2j plays effectively the role of energy as discussed already 
in this work and Part II. Therefore, in view of (5.38), it is appropriate to write: 
fi, = f(Ei), with Ei being the R-G energies. Although the explicit form of such 
/—dependence may be difficult to obtain, for our purposes it is sufficient only 
to know that such a dependence does exist. This then allows us to make an 
identification: Q,i <=^ — - consistent with Varchenko's results, e.g. compare his 

K 

Theorem 3.3.5 (page 35) with Theorem 6.3.2. (page 90) [8]. But, we had estab- 
lished that Cti is an integer, therefore, — should be also an integer. This creates 

K 

some apparent problems. For instance, when |M| = k, the determinant, detM, 
becomes zero implying that solutions of K-Z equation become interdependent. 
This fact has physical significance to be discussed below and in Section 6. To 
do so we use some results from our Part I. In particular, a comparison between 

*fta -f + (1)) ( 5 - 48 ) 



sm7rz = it 2 

fe=i 

and 

— = ze- Cz 

r(z) 

v ; k=l 

where C is some known constant, tells us immediately that not only \M\ = n 
will cause det M =0 but also |M| = n(k + 1), k = 0, 1, 2, ... Accordingly, the nu- 

777- j 

merator of (5.47) will create poles whenever — = 1. Existence of independent 

K 



7 

n( i+ u)) e k ( 5 - 49 ) 

1 — 1 v / 



i ±a= n (^-zj) 

l<i,j<L 
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K-Z solutions is not destroyed if, indeed, such poles do occur. These facts allow 

TJX ' 

us to relabel — - as a(s) (or a(t) or a(u), etc.) as it is done in high energy 
n 

physics with continuous parameters s, t, u,... replacing discrete i's, different 
for different T functions in the numerator of (5.47). In the simplest case, this 
allows us to reduce the determinant in (5.47) to the form used by Lovelace, i.e. 

d etM = -X T 4^^M. (5.50) 
r(l - a{s) - a(t)) v ; 

If, as usual, we parametrize a(s) = a(0) + a's, then equation 1 = a(s) + a(t) 
causes the detM to vanish. This also fixes parameter a(0): a(0) = 1/2. This 
result was obtained by Adler long before sting theory emerged and is known 
as Adler's selfconsistency condition [92]. With such "gauge fixing", one can fix 
the slope a' as well if one notices that the experimental data allow us to make 
a choice: 1 = a(mp).This leads to: a' — ^2 ~ 0.885(Get;~ 2 ) in accord with 
observations.. 

The obtained results arc not limited to study of excitations of just one 
" supeconducting" pair of quarks. In princile, any finite amount of such pairs 
can be studied. In such a case the result for det M becomes considerably more 
complicated but the connections with one dimensional magnets become even 
more explicit. We plan to discuss these issues in future publications. 

6 Discussion. Unimaginable ubiquity of Veneziano- 
like amplitudes in Nature 

6.1 General remarks 



In the Introduction, following Heisenberg, we posed the question: Is combi- 
natorics of observational data sufficient for recovery of the underlying unique 
microscopic model? That is, can we have the complete understanding of such a 
model based on information provided by combinatorics? As we demonstrated, 
especially in Section 4, this task is impossible to accomplish without imposing 
additional constraints which, normally, are not dictated by the combinatorics 
only. In Section 5 we demonstrated that, even accounting for such constraints, 
the obtained results could be in conflict with rigorous mathematics. Last but 
not the least, since Veneziano amplitudes gave birth to string theory one can 
pose a question: Is these Veneziano (or Veneziano-like) amplitudes, perhaps 
corrected to account for particles with spin, contain enough information (an- 
alytical, number-theoretic, combinatorial, etc.) that allows restoration of the 
underlying microscopic model uniquely? The answer is: "No"! In the rest of 
this section we explain why. 
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6.2 Random fragmentation and coagulation processes and 
the Dirichlet distribution 



We begin by recalling some known facts from the probability theory. For in- 
stance, we recall that the stationary Maxwell distribution for velocities of par- 
ticles in the gas is just of Gaussian-type. It can be obtained as a station- 
ary solution of the Boltzmann's dynamical equation maximizing Boltzmann's 
entropjffl-The question arises: Is it possible to find (discrete or continuous) 
dynamical equations which will provide known probability laws as stable sta- 
tionary solutions? This task will involve finding of dynamical equations along 
with the corresponding Boltzmann-like entropies which will reach their max- 
ima at respective equilibria for these dynamical equations. We are certainly not 
in the position in this closing section of our paper to discuss this problem in 
full generality. Instead, we focus our attention only on processes wihich are de- 
scribed by the so called Dirichlet distributions. These originate from the integral 
(equation (2.8) of Part I) attributed to Dirichlet, that is 

V(xi, ...,x n+ i) = J "J ' ' ~ u i u n ) x "+ 1 ~ 1 dui ■ ■ -du n . 

ui>0,..., ti„>0 

?ilH \-U n <l 

(6.1) 

n 

A random vector (Xi, ...,X„) £ R" such that X.; > Vi and X^ =1 is 

i=i 

said to be Dirichlet distributed with parameters (x\, x n ; x n+ i) [nevzorov] if 
the probability density function for (Xi, ...,X„) is given by 

r(zi) • • • v(x n+1 ) ^ 

= ~F1 — ^ rv T u i ■■■< u n+i , provided that 

1 [Xi) ■ ■ ■ 1 (Xn+l) 
= 1 - Ui - ■ ■ ■ - U n 



To get some physical feeling of just defined distribution, we notice the following 
peculiar aspects of this distribution. First, for any discrete distribution, we 
know that the probability pi must be normalized, that is ^2 i Pi = 1- Thus, 
the Dirichlet distribution is dealing with averaging of the probabilities! Or, 
better, is dealing with the problem of effectively selecting the most optimal 
discrete probability. The most primitive of these probabilities is the binomial 
probability given by 

)p m (l-p) n - m , m = 0,1,2,...., n. (6.3) 

As discussed recently in our work [93] on the Poincare' and geometrization conjectures. 
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If X is random variable obeying this law of probability then, the expectation 
E{X) is calculated as 

n 

E(X) = m Pm = np = (i. (6.4) 

m— 1 

Consider such a distribution in the limit: n — > oo. In this limit, if we write 
p = fi/n , then the Poisson distribution is obtained as 

m 

Pm = (6.5) 
to! 

Next, we notice that to! = T(m + 1), furthermore, we replace to by real valued 
variable a and fi by x. This allows us to introduce the gamma distribution with 
exponent a whose probability density is 

Px(x) = Y^)^ 1 ^ (6-6) 

for some gamma distributed random variable X. Finally, we would like to 
demonstrate now how the Dirichlet distribution can be represented through 
gamma distributions. Since the gamma distribution originates from the Pois- 
son distribution, sometimes in literature the Dirichlet distribution is called the 
Poisson-Dirichlct (P-D) distribution [14]. To demonstrate connection between 
the Dirichlet and gamma distributions is relatively easy. Following Kingman 
[14], consider a set of positive independent gamma distributed random vari- 
ables: Yi, Y n+ i with exponents ct\, a n +i- Furthermore, consider Y = 
Y\ + ■ ■ ■ + Y n+ i and construct a vector u with components: Ut = Then, 
since Y^i=i u i = 1j the components of this vector are Dirichlet distributed and, 
in fact, independent of Y. Details are given in Appendix C. 

Such described Dirichlet distribution is an equilibrium measure in various 
fields ranging from spin glasses to computer science, from linguistics to genet- 
ics, from forensic science to economics, etc. Many useful references involving 
these and other applications can be found in [9-11]. Furtheremore, most of 
fragmentation and coagulation processes involve the P-D distribution as their 
equilibrium measure. Some applications of general theory of these processes to 
to nuclear and particle physics were initiated in a series of papers by Mekjian, 
e.g. see [12] and references therein. To avoid duplications, we would like to 
rederive some of Mekjian results differently in order to exibit their connections 
with previous sections. 

6.3 The Ewens sampling formula and Veneziano ampli- 
tudes 

This formula was discussed by Mekjian in [94] without any reference to P-D 
distribution. It is discussed in many other places, including Ewens own mono- 
graph [95]. Our exposition follows work by Watterson [96] where he considers 
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a simple P-D average of monomials of the type generated by the individual 
terms in the expansion^! 

u" = («! + •• • + «*)"= V - n — ^ r< (6-7) 

n=\n\ ,...,nfcj 

Such type of expansion was used in Part I (equations (2. 9), (2. 11)) for calcu- 
lation of multiparticle Veneziano amplitudes. Not surprisingly, Watterson's 
calculation also results in the multiparticle Veneziano amplitude which upon 
multiplication by some combinatorial factor in a well defined limit produces 
the Ewens sampling formula playing a major role in genetics. Although in 
Appendix D we reproduce the Ewens sampling formula without use of the P-D 
distribution, Kingman [98] demonstrated that "A sequence of populations has 
the Evens sampling property if and only if it has the P-D limilr T . Hence, we 
expect that our readers will consult the Appendix D prior to reading of what 
follows. Furthermore, since the vector u is P-D distributed, it is appropriate 
to mention at this point that equation (6.7) genetically represents the Hardy- 
Weinberg law [95] for mating specie^]. Hence, the Ewens sampling formula 
provides a refinement of this law accounting for mutations. 

Considear a special case of (6.2) for which x\ = x-z = • • -xk+i — s and let 
e = 9/ K with parameter 9 to be defined later. Then, (6.2) is converted to 

1?((K -\- l)s) 

Px u ...,x. K (u 1 ,...,u K ) = fe (u) = K+1 Y[ uf" 1 provided that 1 

ll \ £ )\ i=l 

= E • n m * ( 6 - 8 ) 

In view of (6.7), consider an average P(n\, ...,tik) over the simplex A (defined 

by X^t 1 u * = !) § iven b y 



P(ni, ...,tik) = — : — r~ f / " ' / u i J ' ' ' u K K 4>k( u )dui ■ ■ ■ du K - (6.9) 



A 

A straightforward calculation produces: 



n! T((K + l)e) A Tfe + rii) 

F '"—^» -■-„,! [r (g)1 K t V n r((A .' +1)£ ; n) (o.io) 

to be compared with (5.47). Evidently, the parameter e can be identified with k 
in (5.47) and, if we select 9 to be a positive integer, then by replacing n^s with 



33 Very recently Watterson's results were successfully applied to some problems in economics 
[97]. 

34 That is to say, that the Ewens sampling formula implies the P-D distribution and vice 
versa. In the context of high energy physics it is appropriate to mention that the law of con- 
servation of energy-momentum reflected in (2.4) leads to the P-D distribution or, equivalently, 
to the Veneziano formula for multiparticle amplitudes. 

35 E.g. see Wikipedia where it is known as Hardy- Weinberg principle. 
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—n^s we reobtain back (5.47) (up to a constant). To obtain the Ewens sampling 
formula (equation (D.6)) from (6.10) few additional steps are required. These 
are: a) we have to let K — > oo while allowing many of n^s in (6.7) to become zero 
(this explains meaning of the word " sampling" ) , b) we have to order remaining 
n^s in such a way that nm > n(2) > • • * > n (k) > 0,0, ...,0, c) we have to 
cyclically order the remaining n^s in a way explained in the Appendix D by 
introducing c[s as numbers of remaining n'^s which are equal to i. That is we 

have to make a choice between representing r = n (i) or r — S[=i under 

condition that k = Y^i=i c i> °0 finally, just like in the case of Bose (Fermi) 
statistics, we have to multipy the r.h.s. of (6.10) by the obviously looking 
combinatorial factor M = K\/[{c\\ ■ ■ ■ c r l)((K — k)l]. Under such conditions 
wc obtain: T((K + l)e) ~ T{0), T((K + l)e + r) = T(9 + r),^^ = ^. 

Less trivial is the result: Kl/[(K - k)\ [T(e)] k ] -> 6 k . Evidently, the factor 

Til t\ 

— : — : in (6.10) now should be replaced by . Finally, a 

m!n 2 ! • • • n K \ n^y ■ ■ n^) 

moment of thought causes us to replace s by i c %^ in order to arrive at the 

Ewens sampling formula: 

r! r 9 Ci 

P(k;n (1) ,...,n (k) ) = ^ JJ — - (6.11) 

in agreement with (D.6). This derivation was made without any reference to ge- 
netics and is completely model- independent. To demonstrate connections with 
high energy physics in general and with Veneziano amplitudes in particular, we 
would like to explain the rationale behind this formula using absolute minimum 
facts from genetics. 

Genetic information is stored in genes. These are some segments (locuses) 
of the double stranded DNA molecule. This fact allows us to think about the 
DNA molecule as a world line for mesons made of a pair of quarks. Phe- 
nomenologically, the DNA is essentially the chromosome. Humans and many 
other species are diploids. This means that they need for their reproduction 
(meiosis) two sets of chromosomes-one from each parent. Hence, we can think 
of meiosis as process analogous to the meson-meson scattering. We would like 
to depict this process graphically to emphasize the analogy. Before doing so we 
need to make few remarks. First, the life cycle for diploids is rather bizarre. 
Each cell of a grown up organism contains 2 sets of chromosomes. The maiting, 
however, requires this rule to be changed. The gametes (sex cells) from each 
parent carry only one set of chromosomes (that is such cells are haploid !). The 
existence of 2 sets of chromosomes makes individual organism unique because 
of the following. Consider, for instance, a specific trait, e.g. "tall" vs "short". 
Genetically this property in encoded in some genj^l. A particular realization 
of the gene (causing the organism to be, say, tall) is called "allele" . Typically, 
there are 2 alleles -one for each of the chromosomes in the two chromosome set. 

36 This is so because the ci numbers count how many of n'^s are equal to i. 
37 Or in many genes, but we talk about a given gene for the sake of argument. 



44 



For instance, T and t (for "tall" and "short"), or T and T or t and t or, finally, 
t and T (sometimes order matters) . Then, if father donates 50% of T cells and 
50% of t cells and mother is doing to do the same, the offspring is likely going to 
have either TT composition with probability 1/4, or tt (with probability 1/4) 
or tT (with probability 1/4) and, finally, tt with probability 1/4. But, one of 
the alleles is usually dominant (say, T) so that we will see 3/4 of tall people 
in the offspring and 1/4 short. What we just described is the essence of the 
Hardy- Weinberg law based, of course, on the original works by Mendel. Details 
can be found in genetics literature [95]. 

Let us concentrate our attention on a particular locus so that the genetic 
character(trait) of a particular individual is described by specifying its two genes 
at that locus. For N individuals in the population there are 2N chromosomes 
containing such a locus. For each allele, one is interested in knowing the propor- 
tion of 2N chromosiomes at which the gene is realized as this allele. This gives a 
probability distribution over the set of possible alleles which describes a genetic 
make-up of the population (as far as we are only looking at some specific locus). 
The problem now is to model the dynamical process by which this distribution 
changes in time from generation to generation accounting for mutations and 
selection (caused by the environment). Mutation can be caused just by chane 
of one nucleotide along the DNA strancp. Normally, the mutant allele is inde- 
pendent of its parent since once the mutation took place it is very unlikely that 
the corrupt message means anything at all. Hence, the mutant can be either 
"good" (fit) or "bad" (unfit) for life and its contribution can be ignored. If 
u is the probability of mutation per gene per generation then, the parameter 
8 = 4Nu in (6.11). With this information , we are ready to restore the rest of the 
genetic content of Watterson's paper [96]. In particular, random P-D variables 
Xi,X2, ...,X/f denote the allele relative frequences in a population consisting 
of K alleles. Evidently, by construction, they are Dirichlet-distributed. Let 
K — > oo and let k be an experimental sample of representative frequencies 
k <C K. The composition of such a sample will be random, both because of the 
nature of the sampling process and because the population itself is subject to 
random fluctuations. For this reason we averaged the Hardy- Weinberg distribu- 
tion (6.7) over the P-D distribution in order to arrive at the final result (6.11). 
This result is an equilibrium result. Its experimental verification can be found 
in [ewens, watterson2]. It is of interest to arrive at it dynamically along the 
lines discussed in Section 6.2. This is accomplished in the next subsection but in 
a different context. Based on the facts just discussed and comparing them with 
those of Section 2 and Part II, it should be clear that both genetics and physics 
of meson scattering have the same combinatorial origin. All random processes 
involving decompositions r — ^ k =1 Jim (or r — ^ r = i * c «) are the P-D processes 
[9]. 

To conclude this subsection, we would like to illustrate graphically why ge- 
netics and physics of hadrons have many things in common. This is done with 

38 The so called "Single Nucleotide Polymorphism" (SNP) which is detectable either elec- 
trophorctically or by DNA melting experiments, etc. 
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Figure 2: The simplest duality diagram describing meson-meson scattering [99]. 
The same picture describes "collision" of two parental DNA's during meiosis 
and can be seen directly under the electron microscope. E.g. see Fig. 2. 3 in [100], 
page 18. 



help of the figures 2 through 4. 



6.4 Stochastic models for second order chemical reaction 
kinetics involving Veneziano-like amplitudes 

In Section 4 and Appendix A we demonstrated the impotant role of the ASEP 
in elucidating the correct physics. Historically, however, long before the ASEP 
was formulated, the role of stochastic processes in chemical kinetics was already 



* *- 



Figure 3: Non -planar loop Pomeron diagram for meson-meson scattering [101]. 
The same diagram describe homologous DNA recombination, e.g. see fig. 2. 2 in 
[100], page 17. 
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Figure 4: Planar loop meson-baryon scattering duality diagram. The same 
diagram describes the interaction scattering between the triple and double 
stranded DNA helices [102] 

recognized. A nice summary is contained in the paper by McQuarrie [103]. The 
purpose of this subsection is to connect the results in chemical kinetics with 
those in genetics in order to reproduce Veneziano (or Veneziano-like) amplitudes 
as an equilibrium measures for the underlying chemical/biological processes. 

h 

Following Darvey et al [104] we consider a chemical reaction A+B <^ C+D 

k-i 

analogous to the meson-meson scattering processes which triggered the discovery 
of the Veneziano amplitudes. Let the respective concentrations of the reagents 
be a, b, c and d. Then, according to rules of chemical kinetics, we obtain the 
following " equation of motion" 

^ = -kiab + fc_icd. (6.12) 

This equation has to be supplemented with the initial condition. It is obtained 
by accounting for the mass conservation. Specifically, let the initial concentra- 
tions of reagents be respectively asa = ^4(0), (3 — B(0), 7 = (7(0) and S = D(0). 
Then, evidently, a + f3 + 7 + 5 = a + b + c + d, provided that for all times 
a > 0,6 > 0,c > and d > (to be compared with equations (2.1), (2.3)). 
Accounting for these facts, equation (6.12) can be rewritten as 

^ = (fc_ 1 -fc 1 )a 2 -[fc 1 (/3-a) + fc_ 1 (2a + 7 + (5)]a + fc_ 1 (a + 7 )(a + /3). (6.13) 

Theis is thus far is standard result of chemical kinetics. The new element 
emerges when one claims that the the variables a, b, c and d are random but 
are still subject to the mass conservation. Then, as we know already from pre- 
vious subsections, we are dealing with the P-D-type process. New element now 
lies in the fact that this process is dynamical. Following Kingman [105] we 
would like to formulate it in precise mathematical terms. For this purpose, we 
introduce the vector p(t)=(pi(t),..., Pfc(t)) such that it moves randomly on the 
simplex A defined by 

A = {p(t); Pj -> 0,5^^ = 1} (6.14) 
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In our case the possible states of the system at time t which could lead to 
a new state specified by a, 6, c, d at time t + At involving not more than one 
transformation in the time interval At are [104] 

a + l 6+1 c- 1 d- 1 

a- I 6-1 c+1 d+1 ) . (6.15) 
abed 

In writing this matrix, following [104], we assume that random variables o, 6, c 
and d are integers, just like in (2. 3), (2. 4). By analogy with equations of motion 
of Appendix A, using (6.15) we obtain, 

P(a, b, c, d; t + At) - P(a,b,c,d;t) = [fci(a + 1)(6 + l)P(a + 1, b + 1, c - 1, d - l;t) 

+fc_i(c+ l)(rf+l)P(a- 1,6- l,c+l,d+ l;t) 
-{kxab + fc_icd)P(a, 6, c, d; t)]At + (9(At)(6.16) 

In view of the fact that the motion is taking place on the simplex A it is sufficient 
to look at the stochastic dynamics of just one variable, say, a (very much like in 
the deterministic equation (6.13). This replaces (6.16) by the following result: 

j f Pa{t) - fci[(a + l)(a+l + )9-a)P 0+ i(t) + *;_i[(7 + a-a+l)(5 + a-a+l)P _i(t) 

— [k\a{(5 — a + a) + fc_i(7 + a — a) (5 + a — a)]P a (t); provided that 
P„(0) = 1, a = a and P a (0) = if a ^ a. (6.17) 

To solve this equation we introduce the generating function G(x, t) via 

G( X ,t) =Y J Pa{t)x a 

and use this function in (6.17) to obtain the following Fokker-Plank-type equa- 
tion 

^G{x,t) = xil-x^h-xk^-^G+il-x^hiP-a + l) 

+k_ l {2a + 1 + 8- l)x\ — G 
ox 

-k- 1 {a + -i){a + 8){l-x)G{x,t) (6.18) 

This equation admits separation of variables: G(x, t) — S(x)T(t) with solution 
for T(t) in the expected form: T{t) = exp(—\ n kit) leading to the equation for 
S(x) 

x(l-x)(l-Kx)^S(x)^P-a+l+K(2a+^+S-l)x](l-x)^S-[K(a+^){a+p){l-s)-X n ]S(x) = 

(6.19) 

This equation is of Lame- type discussed in Section 5 (e.g. see (5.15)) and, there- 
fore, its solution should be a polynomial in x of degree at most w where w 
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should be equal to the minimum of (a + 7, a + S, (3 + 7, 6 + S). As in quan- 
tum mechanics, this implies that the spectrum of eigenvalues A n is discrete, 
finite and nondegenerate. Among these eigenvalues there must be Ao = since 
such an eigenvalue corresponds to the time-independent solution of (6.19) cor- 
responding to true equilibrium. Hence, for this case we obtain instead of (6.19) 
the following final result: 

d 2 d 
x(l-Kx)—^S(x) + lB-a+l+K(2a+j+S-l)x]— S-[K(a+j)(a+/3)]S = 
dx z dx 

(6.20) 

where K = k-\jk\. This constant can be eliminated from (6.20) if we rescale 
x : x — > Kx. After this, equation acquires the standard hypergeometric form 

d 2 d 
x(l~x) — S(x) + [P-a + l + (2a + j + 6~l)x}—S(x)-(a + ~f)(a + /3)S(x) = 0. 
dx z dx 

(6.21) 

In [105] Kingman obtained the Fokker-Planck type equation analogous to our 
(6.18) describing the dynamical peocess whose stable equilibrium is described by 
(6.21) ( naturally, with different coefficients) and leads to the P-D distribution 
(6.2) essential for obtaining Ewens sampling formula. Instead of reproducing his 
results in this work, we would like to connect them with results of our Section 
5. For this purpose, we begin with the following observation. 

6.4.1 Quantum mechanics, hypergeometric functions and P-D dis- 
tribution 

In our works [2,3] we provided detailed explanation of the fact that all exactly 
solvable 2-body quantum mechanical problems involve different kinds of special 
functions obtainable from the Gauss hypergeometric funcftion whose integral 
representation is given by 

1 

F(a,b,c;z) = T{b ^_ » J " t h -\l-t) c - b -\l- Zt)- a dt. (6.22) 



As it is well known from quantum mechanics, in the case of disctete spec- 
trum all quantum mechanical problems involve orthogonal polynomials. The 
question then arises: under what conditions on coefficients (a, b and c) in- 
finite hypergeometric series whose integral representation is given by (6.22) 
can be reduced to a finite polynomial? This happens, for instance, if we im- 
pose the quantization condition: —a = 0, 1,2, .... In such a case we can write 
(1 — zt)~ a = J27=i(.7 a )(~^y( z ty and use this finite expansion in (6.22). In view 
of (6.2) we obtain the convergent generating function for the Dirichlet distribu- 
tion (6.2). Hence, all known quantum mechanical problems involving discrete 
spectrum are effectively examples of the P-D stochasic processes^^ Next, we 

39 For hypergeometric functions of multiple arguments this was recently shown in [106]. 
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are interested in the following. Given this fact, can we include the determinan- 
tal formula (5.47) into this quantization scheme? Very fortunately, this can be 
done, as explained in the next subsection.. 



6.4.2 Hypergeometric functions, Kummer series expansions and Veneziano 
amplitudes 

In view of just introduced quantization condition, the question arises: is this the 
only condition reducing the hypergeometric function to a polynomial ? More 
broadly: what conditions on coefficients a, b and c should be imposed so that 
the function F(a, b, c; z) becomes a polynomial? The answer to this question 
was provided by Kummer in the first half of 19th century [107]. We would 
like to summarize his results and to connect them with dctcrminantal formula 
(5.43). By doing so we shall reobtain Veneziano amplitudes for chemical process 
described by (6.21). 

According to general theory of hypergeometric equations [107], the infinite 
series for hypergeometric function degenerates to a polynomial if one of the 
numbers 



is an integer. This condition is equivalent to a condition that, at least one of 
eight numbers ±(c— l)±(a— 6)±(a+6— c) is an odd number. According to general 
theory of hypergeometric functions of multiple arguments summarized in Section 
5, the k = 1 -type solutions can be obtained using 1-forms (5.43) accounting 
for singular module constraint (5.45). in the form given by equation (5.42). 
In the case of Gauss- type hypergeometric functions, relations of the type given 
by (5.45) were obtained by Kummer who found 24 interdependent solutions. 
Evidently, this number is determined by the number of independent Pochhamcr 
contours [107]. Therefore, among these he singled out 6 (generating these 24) 
and among these 6 he established that every 3 of them are related to each other 
via equation of the type (5.45). 

Let us denote these 6 functions as ui, ...,ug then, we can represent, say, 112 
and uq using u\ and U5 as basis set. We can do the same with m and u 5 by 
representing them through u 2 and u 6 and, finally, we can connect u 3 and u 4 
with Mi and W5. Hence, it is sufficient to consider, say, ui and uq. We obtain, 



a, b, c — a or c — b 



(6.23) 




(6.24) 
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For instance, we obtain: 

T(a)T(b) _ T(a - c + l)r(6 - c + 1) T(c - a)r(c - 6) _ T(l - a)T(l - b) 

r(c-i) ~ r(i - c ) an r(c-i) ~ r(i - c ) ' 

(6.25) 

For c = 1 we obtain an identity. From [darvey] we find that (6.21) admits 2 
independent solutions: 

, , , either F(-a- 7, -a- 6, 0- a + 1; ifz), for /3 > a , , 

«S(xJ = { 3 (6.26) 

or (if a)""^ F(-/3 - 7, -/3 - £, a - /3 + l\Kx), for /3 < a . 

Hence, the condition c = 1 in (6.25) causes two solutions for S{x) to degenerate 
into one polynomial solution, provided that we make an identification: (3 — a 
in (6.26). Notice that to obtain this result there is no need to impose an extra 
condition: a = J@ (or, in our case, which is the same as 7 = 5). 

This makes sence physically both in chemistry and in high energy physics. 
In the case of high energy physics, if the Veneziano amplitudes are used for 
description of, say, irir scattering, in Part I (page 54) it is demonstrated that 
processes for which " concentrations "a = b cause this amplitude to vanish. The 
Veneziano condition: a + b + c = — 1((1.5) of Part I) has its analog in chemistry 
where it plays the same role, e.g. of mass conservation. In the present case we 
have a + /3 + 7 + <5 = const and the Veneziano-like amplitude obtainable from 
(6. 25), (6. 26) is given now by 

-cT(-c) 

In view of known symmetry of the hypergeometric function: F(a, b, c; x) — 
F(b, a, c; x), we also have: V c (b, a) = V c (a, b). This is compatible with the sym- 
metry for Veneziano amplitude. Combining (5.47) with (6.27) we have the 
following options: a) a — 0,7 = 1,6 = 1,2, ...;&) a = 1,7 = 0,5= 0,1,2,... 
These conditions are compatible with those in (1.19) of Part I for Veneziano 
amplitudes. Finally, in view of (6.22), these are quantization conditions for 
resonances as required.. 



A. Basics of ASEP 



A.l. Equations of motion and spin chains 



The one dimensional asymmetric simple exclusion process (ASEP) had been 
studied for some time [108]. The purpose of this Appendix is to summarize the 
key features of this process which are of immediate relevance to the content of 
this paper. To this purpose, following Schutz [109], we shall briefly describe the 
ASEP with sequential updating. Let Bn := {x±, ...,xn} be a set of sites of one 

40 Here a and b have the same meaning as in (6.22) and should not be confused with con- 
centrations. 
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dimensional lattice arranged at time t in such a way that x\ < x 2 < ■ ■ ■ < x^. 
It is expected that each time update will not destroy this order. 

Consider first the simplest case of N = 1. Let pr (pl) be the probability of 
a particle located at the site x to move to the right (left) then, after transition to 
continuous time, the master equation for the probability P(x; t) can be written 
as follows 

t) = p R P(x - 1; t) + p L P(x + 1; t) - P(x- 1). (A.l) 

Assuming that P(x;t) = cxp(— et)P(x) so that that 

P(x;t) = J dpcxp(—et)f(p)exp(ipx), p <G [0, 2tt), we obtain the dispersion 

o 

relation for the energy e(p) : 

e(p)=p*(l-e-*)+p L (l-e*). (A.2) 

The initial condition P(x; 0) = S XtV determines f(p) = e~ lpy /2ir and yields 
finally 

P{x- 1 i y; 0) = -!- / dpe-^e-^e^ = e- {q+q ^ Dt q x - y I x - y {2Dt), (A.3) 
2tt J 



where q = \/pr/pl, D = ^prPl and I n {2Dt) is the modified Bessel function. 
These results can be easily extended to the case N = 2. Indeed, for this case 
we obtain the following equation of motion 

eP(x 1 ,x 2 ) = -p R (P(x 1 -l,x 2 )+P(x u x 2 -l)-2P(x 1 ,x 2 )) 

-Pl(P(xi + 1, x 2 ) + P{x u x 2 + 1) - 2P( Xl ,x 2 )) (A.4) 

which should be supplemented by the boundary condition 

P(x, x + 1) = p R P{x, x) + p L P(x + 1, x + 1) Vx. (A.5) 

Imposition of this boundary condition allows us to look for a solution of (A.4) 
in the (Bethe ansatz) form 

P(x u x 2 ) = A 12 e ipixl e ip2X2 + A 21 e ip2X1 e ipiX2 (A.6) 

yielding e(pi,p 2 ) = s(pi) +s(p 2 ). Use of the boundary condition (A.5) fixes the 
ratio (the S-matrix) S\ 2 = A\ 2 /A 2 i as follows: 

Q , s p R+PL e* p ^ P2 -e^ 

S(pi,p 2 ) = — — . (A. 7) 

v ' p R +p L e lpi+tp2 - e lp ? y ' 

To connect this result with the quantum spin chains, we consider the case of 
symmetric hopping first. In this case we have pr = pl = 1/2 so that (A. 7) is 
reduced to 

1 _l_ e wi+V2 _ 2e lpi 
Sxxx(pi,P2) = - 1 + eipi+ip2 _ 2eip2 (A.8) 
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from which we can recognize the S matrix for XXX spin 1/2 Hcisenberg fer- 
romagnet [67]. If pr ^ pl, to bring (A. 7) in correspondence with the spin 
chain S— matrix requires additional efforts. Following Gwa and Spohn [110] we 
replace the complex numbers e Wl and e W2 in (A. 7) respectively by z\ and z 2 . 
In such a form (A. 7) exactly coincides with the S-matrix obtained by Gwa 
and 

SpokrS 

After this we can rescale z< (i — 1,2) as follows: z,; = . -Zi. 
Substitution of such an asatz into (A. 7) leads to the result 

SxXz{Zl > Z2) = - l + ~ Zl ~z 2 -2A~z 2 > (A ' 9) 

provided that 2A = 1/ \/plVr. For pr = pl = 1/2 we obtain A = 1 as 
required for the XXX chained If, however, pr pl, then, the obtained S- 
matrix coincides with that known for the XXZ model [goden] if we again relabel 
Zi by e lpi which is always permissible since the parameter p is determined by 
the Bethe equations (to be discussed below) anyway. 

In the case of XXZ spin chain it is customary to think about the massless 
— 1 < A < 1 and massive |A| > 1 regime. The massless regime describes various 
CFT discussed in the text while the massive regime describes massive excitations 
away from criticality. As Gaudin had demonstrated [67], for XXZ chain it is 
sufficient to consider only A > domain which makes XXZ model perfect for 
uses in ASEP. The cases A = and A — > oo also physically interesting: the 
first corresponds to the XY model and the second to the Ising model. 

Once the S-matrix is found, the N— particle solution can be easily con- 
structed [109]. For instance, for N=3 we write 

ty(xi,x 2 ,x 3 ) = exp(ipiXi+ip 2 x 2 +ip3X 3 ) + S 2 iexp(ip 2 xi+ipix 2 +ip3X 3 ) 

+S32S31 ex.p(ip 2 X! + ip 3 x 2 + ipix 3 ) + S 2 iS 31 S3 2 exp(ip 3 xi + ip 2 x 2 + ip 3 x 3 ) 
+S31S32 exp(ip 3 xi + ipix 2 + ip 2 x 3 ) + S 32 exp(ipixi + ip3X 2 + ip 2 x 3 ), 

etc. This result is used instead of f(p) in (A. 3) so that the full solution is given 

by 



N 

P(xi, ...,x N ;t\ y u ...,y N ; 0) - J] — / d Pl e~ £ M t e~^y^(x 1: ...,x N ) 

1=1 11 { 

(A.ll) 

The abobe picture should be refined as follows. First, the particle sitting at 
Xi will move to the right (left) only if the nearby site is not occupied. Hence, 
the probabilities pr and pl can have values ranging from to 1. For instance, 
for the totally asymmetric exclusion process (TASEP) particle can move to the 
right with probability 1 if the neighboring site to its right is empty. Othervice 



41 E.g. see their equation (3.5). 

42 It should be noted though that such a parametrization is not unique. For instance, 
following [56] it is possible to choose a slightly different parametrization, e.g. A = — i (g+g - 1 ) , 

where q = ^/pr/pl- 
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the move is rejected. Since under such circumstances particle can never move 
to the left, there must be a particle source located next to the leftmost particle 
position and the particle sink located immediately after the rightmost position. 
After imposition of emission and absorption rates for these sources and sinks, 
we end up with the Bethe ansatz complicated by the imposed boundary condi- 
tions. Although in the case of solid state physics these conditions are normally 
assumed to be periodic, in the present case, they should be chosen among the 
solutions of the Sklyanin boundary equation [45, 58]. At more intuitive level of 
presentation compatible with results just discussed, the Bethe ansatz for XXZ 
chain accounting for the boundary effects is given in the pedagogically written 
paper by Alcaraz et al [111]. 



A. 2. Dynamics of ASEP and operator algebra 

To make these results useful for the main text, few additional steps are 
needed. For this purpose we shall follow works by Sasamoto and Wadati [112] 
and Stinchcombe and Schiitz [45] . In doing so we rederive many of their results 
differently. 

We begin with observation that the state of one dimensional lattice contain- 
ing N sites can be described in terms of a string of operators D and E, where 
D stands for the occupied and E for empty k— th position along the Id lattice. 
The non normalized probability (of the type given in (A. 11)) can then be pre- 
sented as a sum of terms like this < EEDEDDD ■ ■ ■ E > to be discussed in 
more details below. 

Let C = D + E be the time- independent operator. Then, for the operator 
D to be time-dependent the following commutation relations should hold 

SC + DC = A, (A.12a) 

CS-CD = A, (A.12b) 

DD + DD=[D,S]. (A.12c) 

If A = p L CD -prDC +{p R - p L )D 2 (or A = p L ED - p R DE ), it is possible 
to determine S using equations (A. 12) so that we obtain, 

D=±[A,C- 1 ], (A.13a) 

S=±{A,C~ 1 }, (A.13b) 

provided that 

AC~ 1 D = DC'- 1 A (A.13c) 

with { , } being an anticommutator. 

As before, let us consider the case pr = pl =\ first. This condition leads to 
A = \[C, D]. It is convenient at this stage to introduce an operator D n =C n ~ 1 
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DC n and its Fourier transform T> p —J2 n D n exp(ipn). Using (A. 13a) with A 
just defined leads to the following equation of motion for D„ : 

D-n = \ [Dn+x +D„_i -2D n ] (A.14) 

to be compared with (A.l). Such a comparison produces at once T> p (t) — 
exp(—s(p))"D p (0) so that e(p) = 1 — cosp as before^. Consider now equation 
(A. 13c). Under conditions pn = pl =\ it can be written as 
CDC- l D +DC- X DC =2D 2 or, equivalently, a0 

D n+1 D n + D n £>„-i = 2D n D n . (A.15) 

Following Gaudin [67], we consider a formal expansion D n D m =a exp(ipin + 
ipixn) + (3 exp(ip2^ + ipim) and use it in the previous equation in order to 
obtain: 

aexp(ipi(n + 1) + ip 2 n) + [3 exp(ip 2 {n + 1) + ipm) 
+aexp(ipin + ipi{n — 1)) + j3 exp{ip2ti + ip\(n — 1)) 
= 2aexp(ipin + ip2n) + 2/3 exp(ip2n + ip\n). (A. 16) 

From here we also obtain: 

(aexp(ip 1 n+ip2n)(exp(ipi)+exp(—ip 2 )—2)+f3 exp(ip 1 n+ip2n)(exp(ip 2 )+exp(—ip 1 )—2) — 
and, therefore, 

or \- a 1 + exp(z(pi +p 2 )) ~ 2exp(i Pl ) / A i7\ 

S(p u p 2 ) = - = -— — - — rr — —— exp(t(p2 - P!)) (A.17) 

p 1 + exp(z(pi +P2)j - 2exp{ip 2 ) 

to be compared with (A. 8). An extra factor exp(i(p2 — Pi)) can be actually 
dropped from the S— matrix in view of the following chain of arguments. 
Introduce the correlation function as follows 



-rTVl 



P( Xl ,...,x N ;t i y 1 ,...,y N ;Q)=Z- 1 Tr[D 1 (t)---D N (t)C T 

N 

' dp l e~ e ^ t e~^ y ^(p 1 ,...,p N ), (A.18) 



N 1 

lis 

i=i 



where t>( Pl , ...,p N ) = Z^-Tr[V Px {Q) ■ ■ ■ V PN {t)C N ] and Z N = tr[C N }. In 
arriving at this result the definition of T> p (t), was used along with the fact that 
CDpC" 1 = e~ lp V p . Also, the invariance of the trace under cyclic permuta- 
tions and the translational invariance of the correlation function implying that 
VP(pi, Pn) 7^ only if ^ pi = was taken into account. These conditions 
are sufficient for obtaining the Bethe ansatz equations 

N 

exp(i Pl N) = J] S( PuPj ) V« ^ j, (A.19) 



43 E.g. see (A.2) with p L = p R = 1/2 

44 Since using definition of D n we have: -D n _|_i = CD n C _1 and D n _i = C — 1 D n C. 
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where S(pi,Pj) is the same S— matrix as in (A. 17), except of the missing factor 
exp(i(pi — pj)) which is dropped in view of translational invariancJ^I. 

Extension of these results to the case pr 7^ pi is nontrivial. Because of 
this, we would like to provide some details not shown in the cited references. In 
particular, contary to claims made in [112], we would like to demonstrate that 
the system of equations (A. 12) obtained in [45] is equivalent to the system of 
equations 

[C, S] = 0, (A.20a) 

CD + CT - SD = -p L CD +p R DC, (A.20b) 

DC + DS — TC = plCD — prDC, (A.20c) 

DD + DD = [T, D] (A.20d) 

obtained in [112] with the purpose of describing asymmetric processes. 

To make a comparison between (A. 12) and (A. 20) we notice that (A. 20) has 
operators S and T which cannot be trivially identified with those present in 
(A. 12). Hence, the task lies in making such an identification. For this purpose 
if we assume that S in (A. 20) is the same as in (A. 12) then, in view of (A. 20a), 
by subtracting (A. 12b) from (A. 12a) we obtain: 

DC + CD = 0. (A.21) 

This leads to either D = C^DC or D — -CDC^ 1 . Therefore, taking into ac- 
count that, by construction, C is time-independent, we obtain: D = —CDC~ 1 + 
O, where is some diagonal time- independent matrix operator. 

Next, using these results we multiply (A. 20b) from the right by C -1 and 
(A. 20c) by C _1 from the left, and add them together in order to arrvie at 
equation (19) of [sasam], i.e. 



2D = prC^DC + plCDC- 1 - { PR + p L )D. (A.22) 

Also, by multiplying this result from the right by D we obtain equation (20) of 
[112], that is 

= prDC- x DC + plCDC^D - {pr+ Pl )D 2 , (A.23) 

provided that [T, D] = 0. That this is indeed the case can be seen from the 
same reference where the following result for T is obtained: 

2T= {2+pR-p L )D+p R C- 1 DC -p L CDC-\ (A.24) 

Using it, we obtain: [T, D] = 0, in view of the fact that [C~ 1 DC 1 D] = and 
\CDC~ 1 , D] = since D — —CDC^ 1 + O as we have already demonstrated. 
Furthermore, (A.24) can be straightforwardly obtained by subtracting (A. 20c) 

45 Surely, in case when the effects of boundaries should be accounted, this factor should be 
treated depending on the kind of boundary conditions imposed. 
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(multiplyed by C _1 from the rihgt) from (A. 20a) (multiplycd by C _1 from the 
left). Thus, contary to the claims made in [112], equations (A. 12) and (A. 20) 
are, in fact, equivalent. Nevertheless, as claimed in [112], the system of equations 
(A. 20) is easier to connect with the Bethe ansatz formalism. 

Indeed, using already known fact that D n =C n ~ 1 DC~ n equation (A. 22) 
can be brought into the form: 

Dn = -\p R D n+x + p L D n -i - (pr+pl)D„}. (A.25) 

This result is formally in agreement with previously obtained (A. 14) for the 
fully symmetric case. The authors of [112] have chosen such a normalization 
for probabilities pr and pl that for symmetric case pr = pl = 1 (instead of 
PR. — Pl = 1/2)- To restore the normally accepted condition pr = Pl = 1/2 
requires only to rescale time appropriately. This observation is consitent with 
the fact that the analog of equation (A. 15) (which plays the central role in the 
Bethe ansatz-type calculations) obtained with help of (A. 23) is given by 

PLD n+1 D n +p R D n D n -i = (p R +p L )D n D n (A. 26) 

which holds true wether we choose pr = pl = 1 or p R = p L = 1/50 Obtained 
results allow us to reobtain the S— matrix for the XXZ model in a way already 
described. 

A. 3. Steady- state and q- algebra for the deformed harmonic os- 
cillator 

Using (4.55) we have 

p R DE-p L ED = ((D + E) (A.27) 

Let now D = A\ + B\a and E = A 2 + B 2 a + where Ai and Bi, i = 1,2, are 
some c-numbers. Substituting these expressions back to (A.27) we obtain the 
following set of equations 

C(Ai + At) -eA x A 2 = C, (A.28a) 
where C is some constant to be determined below, and 

CBi = eBxAi, (A.28b) 

(B 2 = eB 2 A 1 , (A.28c) 

From here we obtain: A\ = A 2 = A = C/e, with B\, B 2 being yet arbitrary 
c-numbers and e — pr — pl ■ We can determine these numbers by comparing 

46 It should be noted though that the authors of [112] have erroneously obtained (e.g. see 
their equation (23)) PrD^ +PlD^ +1 = (p R + p L )D n D n+ i instead of our (A. 26). 
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our results with those in [47]. This allows us to select B\ = B 2 = ^yf_ q > 

— = -p— = C, q = — so that we obtain: 
£ 9 PR 

D = -J— + -7=^, (A.29a) 
1 - g V 1 - 9 

£ = , + 7= a + (A.29b) 

1 - g VI - 9 



and, finally, 

in accord with (4.28d). 



aa + - ga+a = 1 (A. 29c) 



B. Linear independence of solutions of K-Z equation 

Linear independence of solutions of K-Z equation is based on the following 
arguments. Consider change of the basis 

&=A>e i , i,j = 1,2,. ..,n (B.l) 

in R™. Using this result, consider the exterior product 

e 1 A • • • Ae" = [dctAJe 1 A • • • Ae" (B.2) 

Next, suppose, that the vectors &> are lineraly -dependent. In particular, this 
means that 

e" = a 1 e 1 + --- + a„_ 1 e"- 1 (B.3) 
for some nonzero a^s. Using this expansion in (B.2) we obtain 

e 1 A • • • A e"" 1 A (c^e 1 + • • • + a„_ie™ _1 ) = (B.4) 

implying [det A] =0. Convesely, if [det A] ^0 then, vectors are linearly in- 
dependent. 



C. Connections between the gamma and Dirichlet distributions 

Using results of our Part I, especially, equation (3.27), such a connection can 
be easily established. Indeed, consider n+ 1 independently distributed random 
gamma variables with exponents a\, a n +i- The joint probability density for 
such variables is given by 

PYl ,... Yn+li su Sn+l) = ^ . . . T^^r- 1 • • • cr - 1 - (CD 
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Let now Sj = Ut, where U are chosen in such a way that En=i*» = L Then > 
using such a substitution into (C.l) wc obtain at once: 

oo 

p„ 1 ,... UB+1 (ti,...,t n +i) = [/ t^e-*]— ■ ■ ■ — 1 t" 1 - 1 ■ ■ ■ Ck" 1 Provided that 1 

J t(aij J- {a n +i) 

o 

x — ^ + 1 

= E !**■ ( c - 2 ) 

* — ^n— 1 

oo 

Since a = ai + ■ ■ ■ + a n +i, we obtain: J i Q_1 e _t = r(«i + • • • + a n+ i) so that 

o 

the density of probability (C.2) is indeed of Dirichlet-type given by (6.2). 
D. Some facts from combinatorics of the symmetric group S n 

Suppose we have a finite set X. \/x e X consider a bijection X — >X made of 

some permutation sequence: x, ir(x), ir 2 {x), ... Because the set is finite, we must 
have n m (x) = x for some m > 1. A sequence (x, n(x),n 2 (x), ..,ir m ~ 1 (x)) = C m 
is called a cycle of length m. The set X can be subdivided into disjoint product 
of cycles so that any permutation it is just a product of these cycles. Normally 
such a product is not uniquely defined. To make it uniquely defined, we have 
to assume that the set X is ordered according to a certain rule. The, standard 
cycle representation can be constructed by requiring that a) each cycle is writ- 
ten with its largest element first, and b) the cycles are written in increasing 
order of their respective largest elements. Let N be some integer and consider a 
decomposition of TV as N — Y^f=i n i ■ We say that n =(n , nx) is partition 
of N (or n h N) The same result can be achieved if, instead we would con- 
sider the following decomposition of N: N — where, according to our 
conventions, we have Ci = c^tt) is the number of cycles of length i. The total 
number of cycles then is given by K=J2iLi c i ■ Define a number S(N, K) as the 
number of permutations of X with exactly K cycles. Then, the Stirling number 
of the first kind can be defined as S{N,K) := {-1) N ~ K S{N, K). The numbers 
S(N, K) can be obtained recursively using the following recurrence relation 

S(N,K) = (N-1)S(N-1,K) + S(N-1,K-1), N,K>1; 5(0,0) = 1. (D.l) 

Use of this recurrence allows us to obtain the following important result 

JV 

§{N, K)x K = x{x + \){x + 2) • • • (x + N - 1). (D.2) 

Let now x = 1 in {D.2), then we can define the probability p(K; N) = S(N, K)/N\ 
Furthermore, one can define yet another probability by introducing a notation 
[a;]^ = x(x + \){x + 2) • • • (x + N - 1). Then, we obtain: 

N K N 

E S(N, K)^ = Pk(N; x) = 1. (D.3) 
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Such defined probability Pk(N; x) can be further rewritten in view of the famous 
result by Cauchy. To obtain his result, we introduce the generating function 

F$ (x) = V ^ -xl 1 ■■■x c " (D.4a) 

and require that S(N, K)x K = JF^(x). This can happen only if 
N W 

^ ^ l c ici!2 c 2 C2 ! • • • N c "c N \ ~~ 1 (D.4b) 

Thus, we obtain 

A" 

/'' <• 



T-r AH ^ w x-^-W 

S(N, K) = [[ -— , provided that K = 2J _ c* and iV = 2J._ zc. (D.5) 



In these notations the Ewens sampling formula acquires the following canonical 
form 



P K (N;x) = X TT iV ", provided that K — a and -/V = ic. 

a; J -- L i Ci cA ^— '*=i ^-^i=i 

L J i= 1 

(D.6) 



A' 



A 



References 



[1] P.Dirac, Lectures on quantum field theory, 

Yeshiva University Press (1996). 
[2] A.Kholodenko, Heisenberg honeycombs solve Veneziano puzzle, 

|hep-th/0608lT7| 

[3] A.Kholodenko, Quantum signatures of Solar System dynamics, 

arXiv.0707.3992. 
[4] A.Kholodenko, New strings for old Veneziano amplitudes I. 

Analytical treatment, J.Geom.Phys.55 (2005) 50. 
[5] A.Kholodenko, New strings for old Veneziano amplitudes II. 

Group-theoretic treatment, J.Geom.Phys.56 (2006)1387. 
[6] A.Kholodenko, New strings for old Veneziano amplitudes III. 

Symplectic treatment, J.Geom.Phys.56 92006) 1433. 
[7] N.Reshetikhin and A.Varchenko, Quasiclassical asymptotics of 
solutions of KZ equations, in Geometry, topology and physics for 
Raoul Bott, p. 293, International Press (1995). 
[8] A.Varchenko, Special functions, KZ type equations, and 
representation theory, AMS Publishers (2003). 



60 



[9] J.Bertoin, Random fragmentation and coagulation processes, 

Cambridge University Press, Cambridge U.K. (2006). 
[10] J. Pitman, Combinatorial stochasic processes, 

Springer- Verlag, Berlin (2006). 
[11] R. Arratia, A.Barbour and S.Tavare, Logarithmic combinatorial 

structures: a probabilistic approach, European Mathematical Society, 

Zurich (2003). 

[12] A.Mekjian, Model for studing brancing processes, multiplicity 

distribution and non-P oissonian fluctuations in heavy -ion collisions, 
PRL 86 (2001) 220. 

[13] R.Stanley, Combinatorics and commutative algebra, Birkhauser, 
Boston (1996). 

[14] J. Kingman, Poisson processes, Clarendon Press, Oxford (1993). 
[15] S.Ghorpade and G.Lachaud, Hyperplane sections of Grassmannians 
and the number of MDS linear codes, Finite Fields & 
Their Applications 7 (2001) 468. 
[16] R.Stanley, Enumerative combinatorics, Vol.1, 

Cambridge University Press, Cambridge, U.K. (1999). 
[17] S.Mohanty, Lattice path counting and applications, 

Academic Press, New York (1979). 
[18] R.Bott and L.Tu, Differential forms in algebraic topology, 

Springer- Verlag, Berlin (1982). 
[19] J.Schwartz, Differential geometry and topology, Gordon 

and Breach, Inc., New York (1968). 
[20] M.Stone, Supersymmetry and quantum mechanics of spin, 

Nucl.Phys. B 314 (1989) 557. 
[21] O.Alvarez, I Singer and P.Windey, Quantum mechanics and 
the geometry of the Weyl character formula, 
Nucl.Phys.B337 (1990) 467. 
[22] A.Polyakov, Gauge fields and strings, 

Harwood Academic Publ., New York (1987). 
[23]. A.Polychronakos, Exact spectrum of SU(n) spin chain with 

inverse square exchange, Nucl.Phys. B419 (1994) 553. 
[24] H. Frahm, Spectrum of a spin chain with invese square exchange, 

J.Phys.A 26 (1993) L473. 
[25] A.Polychronakos, Generalized statistics in one dimension, 
|hep-th/9902157| 

[26] A.Polychronakos, Physics and mathematics of Calogero particles, 
|hep-th/0607033| 

[27] K.Hikami, Yangian symmetry and Virasoro character in a lattice 

spin system with long -range interactions, 

Nucl.Phys.B 441 (1995) 530. 
[28] E.Melzer, The many faces of a character, hep-th/9312043. 
[29] R.Kedem, B. McCoy and E.Melzer, The sums of Rogers, Schur and 

Ramanujian and the Bose-Fermi correspondemnce in 1+1 

dimensional quantum field theory, hep-th/930^056 . 



61 



[30] A.Tsvelik, Quantum field theory in condensed matter physics, 

Cambridge University Press, Cambridge U.K. (2003). 
[31] J.Goldman and J. C. Rota, The number of subspaces of a vector space, 

in Recent Progress in Conbinatorics, 

p. 75, Academic Press, New York (1969). 
[32] V.Kac and P.Cheung, Quantum calculus, Springer- Verlag, Berlin (2002). 
[33] G.Andrews, The theory of partitions, 

Addison- Wesley Publ.Co., London (1976). 
[34] D.Galetti, Realization of the q-deformed harmonic oscillator: 

Rogers-Szego and Stiltjes- Wiegert polynomials, 

Brazilian Journal of Physics 33 (2003)148. 
[35] M.Chaiichian, H.Grosse and P.Presnajer, Unitary representations 

of the q-oscillator algebra, J.Phys.A 27 (1994) 2045. 
[36] R. Floreani and L.Vinet, Q-orthogonal polynomials and the oscillator 

quantum group, Lett.Math.Phys. 22 (1991) 45. 
[37] H.Karabulut, Distributed Gaussian polynomials as q-oscillator 
eigenf unctions, 

J.Math.Phys.47 (2006) 013508. 
[38] A.Macfarlane, On q-analogues of the quantum harmonic oscillator 

and the quantum group SU(2) q , J.Phys.A 22 (1989) 4581. 
[39] H.Karabulut and E.Siebert, Distributed Gaussian polynomials 

and associated Gaussian quadratures, J.Math.Phys.38 (1997) 4815. 
[40] A.Atakishiev and Sh.Nagiyev, On the Rogers-Szego polynomials, 

J.Phys.A 27 (1994) L611. 
[41] G. Gasper and M.Rahman, Basic hypergeometric series, 

Cambridge University Press, Cambridge, U.K. (1990). 
[42] R.Koekoek and R.Swarttouw, The Askey- scheme of hypergeometric 

orthogonal polynomials and its q-analogs, arXiv:math/9602214 
[43] M. Ismail, Classical and quantum orthogonal polynomials of one 

variable, Cambridge University Press, Cambridge, U.K. (2005). 
[44] T.Nagao and T.Sasamoto, Asymmetric simple exclusion process and 

modified random matrix ensembles, Nucl.Phys.B 699 (2004) 487. 
[45] R.Stinchcombe and G.Schutz, Application of operator algebras to 

stochastic dynamics and Heisenberg chain, PRL 75 (1995) 140. 
[46] T.Sasamoto, One- dimensional partially asymmetric simple exclusion 

process with open boundaries: orthogonal polynomials approach, 

J.Phys.A 32 (1999) 7109. 
[47] R.Blythe, M.Ewans, F.Colaiori and F.Essler, Exact solution of 

a partially asymmetric exclusion model using a deformed 

oscillator algebra, J.Phys.A 33 (2000) 2313. 
[48] B.Derrida, M.Ewans, V. Hakim and V.Pasquer, Exact solution 

of a Id asymmetric exclusion model using a matrix formulation, 

J.Phys A 26 (1993) 1493. 
[49] B.Derrida and K.Malick, Exact diffusion constant for the one- 
dimensional partially asymmetric exclusion model, 

J.Phys.A 30 (1997) 1031. 



62 



[50] M.Kardar, G.Parisi and Yi-Ch. Zhang, Dynamic scaling of groving 

interfaces, PRL 56 (1986) 889. 
[51] T.Sasamoto, S.Mori and M.Wadati, One- dimensional asymmetric 
exclusion model with open boundaries, 

J.Phys.Soc.Japan 65 (1996) 2000. 
[52] T. Oliviera, K.Dechoum, J.Redinz and F. Aarao Reis, Universal 

and nonuniversal features of the crossover from linear to 

nonlinear interface growth, 

Phys.Rev. E 74 (2006) 011604. 
[53] A. Lazarides, Coarse- graining a restricted solid-on -solid model, 

Phys.Rev.E 73 (2006) 041605. 
[54] D.Huse, Exact exponents for infinitely many new multicritical points, 

Phys.Rev.B 30 (1984) 3908. 
[55] D.Friedan, Z.Qui and S.Shenker, Conformal invariance, unitarity 

and critical exponents in two dimensions, PRL 52 (1984) 1575. 
[56] J.de Gier and F.Essler, Exact spectral gaps of the asymmetric 

exclusion process with open boundaries, J.Sat.Mech.(2006) P12011. 
[57] V.Pasquer and H.Saleur, Common structures beetween finite 

systems and conformal field theories through quantum groups, 

Nucl.Phys. B 330 (1990) 523. 
[58] C.Gomez, m.Ruiz-Altaba and G.Sierra, 

Quantum groups in two-dimensional physics, 

Cambridge University Press, Cambridge, U.K. (1996). 
[59] L.Faddeev and O.Tirkkonen, Connections of the Liouville model 

and XXZ spin chain, Nucl.Phys.B453 (1995) 647. 
[60] A.Kholodcnko, Kontsevich-Witten model from 2+1 gravity: 
New exact combinatorial solution, J.Geom.Phys.43 (2002) 45. 
[61] P.Forrester, Visious random walkers in the limit of a large 

number of walkers, J.Stat.Phys. 56 (1989) 767. 
[62] T.Sasamoto, Fluctuations of the one- dimensional asymmetric 

exclusion process using random matrix techniques, 
J.Stat.Mech. (2007) P07007. 
[63] A.Mukherjee and S.Mukhi, c=l matrix models: equivalences 

and open-closed string duality, JHEP 0510 (2005) 099. 
[64] J.Distler and C.Vafa, A critical matrix model at c=l, 

Mod.Phys.Lett.A 6 (1991) 259. 
[65] D.Ghoshal and C.Vafa, c=l strting as the topological theory 

on a conifold, Nucl.Phys. B 453 (1995) 121. 
[66] D.Huse and M.Fisher, Commensurate melting, domain walls, 

and dislocations, Phys.Rev.B29 (1984) 239. 
[67] M.Gaudin, La function d'onde de Bethe, 

Masson, Paris (1983). 
[68] D.Grabiner, Brownian motion in a Weyl chamber, 
non -colliding particles, and random matrices, 
arXiv: |math.RT/9708207| 
[69] C. Krattenhaler, Asymptotics for random walks in alcoves 



63 



of affine Weyl groups, arXiv: math/0301203. 
[70] S.de Haro, Chern-Simons theory, 2d Yang -Mills, 

and Lie algebra wanderers, Nucl.Phys. B730 (2005) 313. 
[71] M.Mehta, Random matrices, Elsevier, Amsterdam (2004). 
[72] J.Maldacena, G.Moore, N.Seiberg and D.Shih, 

Exact vs. semiclassical target space of the minimal string, 
JHEP 0410 (2004) 020. 
[73] K.Okuyama, D-brane amplitudes in topological string on conifold, 

Phys.Lett.B 645 (2007) 273. 
[74] M.Tierz, Soft matrix models and Chern-Simons partition functions, 

Mod.Phys. A 19 (2004) 1365. 
[75] P.Etingof, I.Frcnkel and A.Kirillov Jr., Lectures on representation 
theory and Knizhnik-Zamolodchikov equations, 
AMS Publishers, Providence, R.I. (1998). 
[76] A.Chervov and D.Talalaev, Quantum spectral curves, quantum 
integrable systems and the geometric Langlands correspondence, 
arXiv: |hep-th/0604128| 
[77] E.Frenkel, Langlands correspondence for loop groups, 

Cambridge University Press, Cambridge, U.K. (2007). 
[78] E.Frenkel and E.Witten, Geometric endoscopy and mirror symmetry, 
arXiv: 0710.5939. 

[79] R.Richardson and N.Sherman, Exact eigenvalues of the pairing-force 

Hamiltonian, Nucl.Phys.52 (1964) 221. 
[80] R.Richardson, Exactly solvable many-boson model, JMP 9 (1968) 1327. 
[81] N.Vilenkin, Special functions and theory of group representations, 

Nauka, Moscow (1991). 
[82] J.Dukelsky, S.Pittel and G.Sierra, Exactly solvable 

Richardson- Gaudin models for many-body 

quantum systems, Rev. Mod.Phys. 76 (2004) 643. 
[83] A.Gorsky, Gauge theories as string teories: the first results, 

arXiv: |hep-th/0602184| 
[84] A.Balantekin, T.Dereli and Y.Pehlivan, Exactly solvable pairing model 

using an extension of Richardson- Gaudin approach, 

Int.J.Mod.Phys.E 14 (2005) 47. 
[85] A.Ushveridze, Quasi-exactly solvable models in quantum mechanics, 

IOP Publishing Ltd., Philadelphia (1994). 
[86] A.Ovchinnikov, Exactly solvable discrete BCS-type Hamiltonians 

and the six-vertex model, Nucl.Phys. B707 (2002) 362. 
[87] M.Alford, A.Schitt, K.Rajagopal and Th.Schafer, 

Color superconductivity in quark matter, 

arXiv: 0709.4635. 
[88] L.Cooper, Boubd electron pairs in a degenerate electron gas, 

Phys.Rev.104 (1956) 1189. 
[89] A.Balantekin, J.de Jesus and Y.Pehlivan, Spectra and symmetry in 

nuclear pairing, Phys.Rev. C 75 (2007) 064304. 
[90] C.Lovelace, A novel application of Regge trajectories, 



64 



Phys.Lett.B 28 (1968) 264. 
[91] A.Kholodcnko, Traces of mirror symmetry in Nature, 

International Math.Forum 3 (2008) 151. 
[92] S.Adler, Consistency conditions on a strong interactions 

implied by a partially conserved axial- vector current, 
Phys.Rev. 137 (1965) B1022. 
[93] A.Kholodcnko, Towards physically motivated proofs 
of the Poincare and geometrization conjectures, 
J.Gcom.Phys.58 (2008) 259. 
[94] A.Mckjian, Cluster distribution in physics and genetic 

diversity, Phys.Rcv.A 44(1991) 8361. 
[95] W.Ewens, Mathematical population genetics, 

Springer- Verlag, Berlin (2004). 
[96] G.Watterson, The stationary distribution of the 

infinitely-many neutral alleles diffusion model, 

J.Appl.Probability 13 (1976) 639. 
[97] M.Aoki, Open models of share markets with two 

dominant types of participants, J. of Economic 

Behaviour & Organization 49 (2002) 199. 
[98] J.Kingman, The population structure associated 

with the Ewens sampling formula, 

Theoretical Pop.Biology 11 (1977) 274. 
[99] H.Harari, Duality diagrams, PRL 22 (1969) 562. 
[100] D. Leach, Genetic recombination, 

Blackwell Science Ltd., Oxford, U.K. (1996). 
[101] P.Freund, Two component duality and strings, 

arXiv: 0708.1983. 
[102] S.Mirkin, Structure and biology of H DNA, 

in Triple Helix Forming Oligonucleotides, pl95, 

Kluwer Academic, Boston (1999). 
[103] D.McQuarrie, Stochastic approach to chemical kinetics, 

J.Appl.Prob. 4 (1967) 413. 
[104] I.Darvey, B.Ninham and P.Staff, Stochastic models for 

second-order chemical reaction kinetics. 

The equilibrium state, J.Chem.Phys.45 (1966) 2145. 
[105] J.Kingman, The dynamics of neutral mutation, 

Proc.Roy.Soc. London A 363 (1978) 135. 
[106] A. Lijoi and E.Regazzini, Means of a Dirichlet process 
and multiple hypergeometric functions, 

Ann.Probability 32 (2004) 1469. 
[107] H.Batcman and A.Erdelyi, Higher transcendental functions, 

Vol.1. McGraw Hill, New York (1953). 
[108] R. Blythe and M.Evans, Nonequilibrium steady states of 

matrix-product form: a solver guide, J.Phys.A 40 (2007) R333. 
[109] G.Schiitz, Exact solution of the master equation for the 

asymmetric exclusion process, J.Stat.Phys. 88 (1997) 427. 



65 



L.Gwa and H.Spohn, Bethe solution for the dynamical- scaling 
exponent of the noisy Burgers equation, Phys.Rcv.A 46 (1992) 844. 
F.Alcaraz, M.Barber, M.Batchelor, R. Baxter and G.Quispei, 
Surface exponents of the quantum XXZ, Ashkin-Teller and 
Potts models, J.Phys.A 20 (1987) 6397. 

T.Sasamoto and M.Wadati, Dynamic matrix product ansatz and 
Bethe ansatz equation for asymmetric exclusion process 
with periodic boundary, J.Phys.Soc.Jpn. 66 (1997) 279. 



66 



